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Deriving  Sea  Surface  Salinity  and  Density  Variations 
From  Satellite  and  Aircraft  Microwave  Radiometer 
Measurements:  Application  to  Coastal 
Plumes  Using  STARRS 

Derek  Burrage,  Joel  Wesson,  and  Jerry  Miller 


Abstract — Using  brightness  temperature  Tb  measurements 
from  L-band  airborne  microwave  radiometers,  with  independent 
sea  surface  temperature  (SST)  observations,  sea  surface  salinity 
(SSS)  can  be  remotely  determined  with  errors  of  about  1  psu  in 
temperate  regions.  Nonlinearities  in  the  relationship  between  T6, 

SSS,  and  SST  produce  variations  in  the  sensitivity  of  salinity  S 
to  variations  in  Tb  and  SST.  Despite  significant  efforts  devoted  to 
SSS  remote  sensing  retrieval  algorithms,  little  consideration  has 
been  given  to  deriving  density  D  from  remotely  sensed  SSS  and 

SST.  Density  is  related  to  5  and  T  through  the  equation  of  state. 
It  affects  the  ocean’s  static  stability  and  its  dynamical  response  to 
forcings.  By  chaining  together  two  empirical  relationships  (flat-sea 
emissivity  and  equation  of  state)  to  form  an  inversion  algorithm 
for  sea  surface  density  (SSD)  in  terms  of  Tb  and  SST,  we  de¬ 
velop  a  simple  L-band  SSD  retrieval  algorithm.  We  use  this  to 
investigate  the  sensitivity  of  SSD  retrievals  to  observed  Tb  and 
SST  and  infer  errors  in  D  for  typical  sampling  configurations 
of  the  airborne  Salinity,  Temperature,  And  Roughness  Remote 
Scanner  (STARRS)  and  satellite-borne  Soil  Moisture  and  Ocean 
Salinity  (SMOS)  and  Aquarius  radiometers.  We  then  derive  D 
from  observations  of  river  plumes  obtained  using  STARRS  and 
demonstrate  several  oceanographic  applications:  the  observations 
are  used  to  study  variations  in  T  and  S  effects  on  D  in  the 
Mississippi  plume,  and  the  across-shelf  density  gradient  is  used 
to  infer  surface  geostrophic  shear  and  subsurface  geostrophic 
current  in  the  Plata  plume.  Future  basin-scale  applications  of 
SSD  retrievals  from  satellite-home  microwave  radiometers  such 
as  SMOS  and  Aquarius  are  anticipated. 

Index  Terms — Aquarius,  brightness  temperature,  coastal 
plumes,  density,  emissivity,  microwave  radiometry,  retrieval 
algorithm,  Salinity,  Temperature,  And  Roughness  Remote 
Scanner  (STARRS),  sea  surface  salinity  (SSS),  Soil  Moisture  and 
Ocean  Salinity  (SMOS). 
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I.  Introduction 

RECENT  advances  in  microwave  remote  sensing  enable 
precise  measurements  of  the  brightness  temperature  (Tb) 
of  seawater  at  L-band  (~1.4  GHz)  to  be  obtained  over  coastal 
regions  from  a  light  aircraft  [1],  [2],  whereas  the  Soil  Moisture 
and  Ocean  Salinity  (SMOS)  and  Aquarius  satellite  missions 
to  be  launched  in  the  near  future  will  provide  a  similar  capa¬ 
bility  over  the  deep  global  oceans  [3]— [6].  Modem  airborne 
L-band  microwave  radiometers  are  sensitive  enough  to  detect 
Tb  variations  at  better  than  0.5-K  precision  under  operational 
conditions,  at  typical  spatial  resolutions  of  1  km,  whereas 
satellite-borne  systems  are  being  designed  to  provide  effective 
0.1 -K  precision,  after  data  averaging,  on  spatial  scales  of  about 
100-300  km. 

The  brightness  temperature  at  the  surface  Tb  is  proportional 
to  the  sea  surface  temperature  (SST)  and  the  emissivity.  The 
emissivity  is  itself  a  function  of  SST  and  conductivity  (which 
together  determine  the  salinity),  roughness  (primarily  deter¬ 
mined  by  wind  speed),  beam  incidence  angle  0 ,  and  antenna 
polarization  p.  The  relationship  between  Tb  and  sea  surface 
salinity  (SSS),  SST,  0 ,  and  p  for  a  flat  calm  sea  may  be  ob¬ 
tained  from  the  so-called  “flat  sea”  empirical  emissivity  model 
relationships,  e.g.,  [7]  and  [8],  and  various  semianalytical  or 
empirical  relationships  to  account  for  the  roughness  influence 
are  available  [9],  [10].  By  taking  into  account  reflected  sources 
of  radiation,  both  atmospheric  and  extraterrestrial  (galactic, 
cosmic  background,  and  solar),  as  well  as  atmospheric  absorp¬ 
tion  effects,  it  is  possible  to  edit  or  correct  the  Tb  values  and, 
together  with  SST,  invert  them  to  estimate  SSS.  Based  on  the 
Klein  and  Swift  relation  [7],  a  precision  of  0.5  K  corresponds 
to  a  retrieved  SSS  precision  of  1  psu  (for  single-point  measure¬ 
ments  under  temperature  conditions  typically  encountered  in 
mid-  and  low-latitude  coastal  seas). 

Significant  effort  has  been  expended  to  date  to  define  and 
refine  the  algorithms  for  SSS  retrieval  from  L-band  Tb ,  taking 
into  account  the  aforementioned  factors.  However,  the  possi¬ 
bility  of  deriving  sea  surface  density  (SSD)  indirectly  from 
retrieved  SSS  and  SST  data  or  directly  from  the  SST  and  Tb 
remote  sensing  measurements  has  received  little  consideration. 
This  is  so,  despite  the  fundamental  importance  of  density  in 
governing  ocean  wave  and  circulation  dynamics,  including 
baroclinic  currents  and  shear  instabilities,  convective  overturn. 
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and  various  air-sea  interaction  processes.  Both  the  seawater 
equation  of  state  (SES),  which  determines  D  as  a  function  of 
S ,  7\  and  pressure  P  [1 1  ],  and  the  equations  governing  flat-sea 
emissivity  [7],  which  determine  Tb  from  SST,  SSS,  6 ,  and  p 
are,  in  general,  nonlinear.  As  demonstrated  here,  the  sensitivity 
of  SST,  SSS,  and  SSD  to  Tb  significantly  changes  throughout 
the  range  of  brightness  temperatures  typically  encountered  over 
the  coastal  ocean.  In  the  text,  we  use  the  acronyms  SST, 
SSS,  and  SSD  to  denote  the  sea  surface  temperature,  salinity, 
and  density,  respectively.  The  corresponding  mathematical  or 
graphic  symbols  are  T  (or  Ts,  to  emphasize  the  surface  value 
and  distinguish  it  from  brightness  temperature  Tb ),  S,  and 
D,  where  surface  values  are  understood,  except  where  T,  S, 
or  D  are  implicit  or  explicit  functions  of  depth,  as  in  T(z), 
S(z ),  or  D(z).  Similarly,  P(z)  denotes  the  pressure  at  depth 
2  in  the  ocean,  unless  2  =  0,  in  which  case  it  refers  to  sur¬ 
face  (atmospheric)  pressure.  The  radiometers  actually  measure 
brightness  temperature  Ta  at  the  antenna,  but  radiometer  data 
processing  algorithms  transform  this  to  its  surface  value  Tb  by 
applying  various  environmental  corrections  based  on  a  radiative 
transfer  model. 

We  formulate  a  retrieval  algorithm  for  SSD  using  the  avail¬ 
able  empirical  functions  for  the  equation  of  state  and  flat-sea 
emissivity.  Accounting  for  the  controlling  influence  of  antenna 
incidence  angle  and  polarization,  we  then  compute  the  sensi¬ 
tivity  of  the  derived  SSD  to  SST  and  SSS.  This  allows  the 
SSD  error  for  particular  SST  and  Tb  input  error  levels  to  be 
determined  via  the  emissivity  model  and  the  SES. 

At  first  order,  observed  Tb  varies  directly  with  respect  to 
SST  and  inversely  with  respect  to  SSS,  whereas  SSD  varies 
inversely  with  respect  to  SST  and  directly  with  respect  to  SSS. 
This  implies  an  inverse  relation  between  Tb  and  D  that  is 
simultaneously  determined  by  S  and  Ts.  These  relationships 
may  affect  the  detection  of  surface  density  fronts.  Considering 
a  specific  example,  compensating  effects  on  D  of  variations 
in  S  and  T  are  sometimes  observed  across  upward-sloping 
coastal  fronts  where  warm  saline  ocean  water  outcrops  above 
and  seaward  of  cooler  and  fresher  coastal  water  [12].  In  this 
density-compensating  situation,  both  SSD  and  Tb  variations 
across  the  outcropping  front  will  be  small,  and  the  front  will 
be  better  located  using  SST  and  SSS  retrievals.  In  the  opposite 
situation,  where  S  and  T  reinforce  each  other  in  their  effect  on 
D,  both  Tb  and  SSD  may  vary  strongly,  but  inversely,  across 
the  front,  which  consequently  may  be  more  reliably  detected 
using  SSD  retrieval.  In  general,  the  errors  associated  with  the 
retrieval  of  SSD  and  SSS  from  Tb  may  be  compensating  or 
reinforcing  to  a  lesser  or  greater  extent,  depending  on  the  water 
mass  locations  in  the  T-S  parameter  space.  As  demonstrated 
here,  a  simple  vector  index,  based  on  the  T-  and  S'-dependent 
density  ratio  [13],  can  be  employed  to  determine  which  of 
SSS  and  SST  dominates  the  SSD  variations  and  whether  the 
resulting  T  and  S  fields  are  significant  in  either  compensating 
or  reinforcing  the  density  variation. 

We  have  applied  the  algorithms  for  computing  SSD  to 
measurements  from  the  Salinity,  Temperature,  And  Roughness 
Remote  Scanner  (STARRS)  of  the  Naval  Research  Laboratory 
(NRL)  (see  Section  V  and  Appendix  A  for  a  detailed  de¬ 
scription  of  the  instrument).  The  measurements  were  obtained 


during  separate  surveys  of  the  Mississippi  River  plume  (in  the 
northern  Gulf  of  Mexico)  and  the  La  Plata  River  plume  (off  the 
South  American  Atlantic  Coast).  Both  plumes  were  surveyed  in 
summer  and  winter.  A  variety  of  conditions  illustrate  how  the 
sensitivity  of  SSS  and  SSD  retrievals  to  observational  errors 
varies  and  reveal  the  potential  benefits  of  deriving  SSD  from 
microwave  L-band  data  in  remote  sensing  studies  of  coastal 
river  plumes. 

This  paper  has  several  innovative  features.  To  our  knowledge, 
this  is  the  first  time  that  SSD  has  been  derived  from  remote 
L-band  microwave  radiometer  observations.  Similarly,  the  ex¬ 
ploration  of  the  combined  parameter  space  of  the  equation  of 
state  and  the  L-band  emissivity  model  and  the  presentation 
of  graphs  illustrating  the  sensitivity  of  SSD  to  variations  and 
errors  in  SST,  SSS,  and  Tb  are  novel.  Furthermore,  surface 
geostrophic  currents  have  previously  been  determined  directly 
from  radar  altimeter  measurements  of  sea  surface  height  (SSH) 
by  various  investigators  (e.g.,  [14]  and  [15],  based  on  the  geo- 
strophy  theory;  and  [16],  combining  geostrophic  with  Ekman 
theory  using  auxiliary  wind  data)  and  indirectly  from  thermal 
infrared  (IR)  and  microwave  SST  measurements  (e.g.,  [17], 
based  on  the  quasi-geostrophic  (QG)  theory).  However,  the 
derivation  of  near-surface  geostrophic  shear  and  subsurface 
geostrophic  currents  from  a  combination  of  remotely  sensed 
SSS  and  SST  data  and  thermal  wind  equations  (through  ap¬ 
propriate  assumptions),  has  not  previously  been  demonstrated. 
Finally,  a  simple  extension  of  the  Turner  angle  [  1 8],  [  1 9]  is  pro¬ 
posed  to  give  a  vector  representation  of  both  horizontal  density 
compensation  and  the  combined  strength  of  the  associated  T 
and  S  gradients. 

In  the  remainder  of  this  paper,  we  examine  the  relation¬ 
ships  determining  SSD  from  Tb  and  SST,  develop  the  SSD 
algorithm,  and  investigate  SSD  sensitivity  to  SST  and  Tb 
errors  (Section  II).  We  then  consider  how  density  varies  and 
introduce  the  density  compensation  index,  i.e.,  the  “Turner 
vector”  (Section  III).  The  results  are  used  to  study  horizontal 
surface  density  gradients  and  to  develop  a  method  to  infer 
near-surface  currents  (Section  IV).  We  proceed  to  describe 
relevant  field  results  and  demonstrate  the  application  of  these 
techniques  (Section  V).  Finally,  we  discuss  the  algorithm  and 
possible  enhancements,  evaluate  the  results  (Section  VI),  and 
draw  conclusions  (Section  VII).  Details  of  the  field  survey 
instrumentation  and  logistics  and  mathematical  derivations  are 
provided  in  the  Appendixes. 

II.  Deriving  SSD 
A.  Density  Retrieval  Algorithm 

The  empirical  Klein  and  Swift  algorithm  (KS77)  [7],  based 
on  the  Debye  model  of  the  dielectric  constant  [20],  is  com¬ 
monly  used  to  determine  the  brightness  temperature  Tb  of  a 
calm  (“flat”)  sea  surface,  given  its  physical  temperature  Ts 
and  salinity  S  [Fig.  1(a)].  The  brightness  temperatures  (vertical 
axis)  are  for  vertically  polarized  (V-Pol)  signals,  at  a  represen¬ 
tative  incidence  angle  of  22°,  that  were  chosen  to  coincide  with 
the  intermediate  right  and  left  beams  of  the  STARRS  L-band 
radiometer  that  we  used  for  algorithm  validation  (see  Section  V 
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(a)  (b) 

Density  =(1000  5:1030]  kg  m3,  Salinity=(0:5:40]  psu  Density=(1 000:5: 1030]  kg  nrv3.Bright  Temp=(94  4  122]  K 


Fig.  1 .  Nonlinear  empirical  functions  determining  salinity  S  and  density  D  from  L-band  brightness  temperature  Tb  and  physical  surface  temperature  Ts  at  an 
incidence  angle  of  22°.  (a)  Flat-sea  emissivity  (KS77)  showing  S  (red  or  gray,  dashed)  and  D  (blue  or  black,  solid)  versus  Tb  and  Ts.  (b)  Equation  of  state 
showing  D  (blue  or  black,  solid)  and  Tb  (red  or  gray,  dashed)  versus  Ts  and  S.  Ranges  and  corresponding  increments,  in  the  format  (low:incr:high]  for  S ,  D , 
and  Tb.  appear  at  the  top,  and  the  low  and  high  contour  lines  are  labeled.  The  dotted  lines  trace  relevant  turning  points. 


and  Appendix  A-l  for  instrument  specifications).  For  the  V-Pol 
signals,  the  brightness  temperatures  for  a  given  Ts  and  S  are 
progressively  higher  at  higher  incidence  angles,  but  the  choice 
of  incidence  angle  does  not  otherwise  change  the  qualitative 
features  mentioned  here.  For  horizontally  polarized  (H-Pol)  sig¬ 
nals,  the  brightness  temperatures  are  the  same  as  those  for  V-Pol 
at  nadir  but  are  progressively  lower  at  higher  incidence  angles. 
The  pattern  of  variation  of  5  for  H-Pol  is  otherwise  quite  similar 
to  that  for  V-Pol.  The  KS77  algorithm,  as  originally  defined,  has 
a  validity  range  of  4-35  psu  in  S, 0  °C  to  30  °C  in  T,  and  0°  to 
90°  in  0 ,  and  it  applies  to  radiation  at  frequencies  lower  than  the 
X-band  (including,  in  particular,  the  L-  and  C-band,  which  are 
detected  by  our  microwave  instrumentation).  Since  our  experi¬ 
ence  shows  that  it  produces  plausible  results  for  S  down  to  0 
and  up  to  at  least  40  psu,  we  have  chosen  to  plot  values  for  this 
extended  range.  When  used  as  a  salinity  retrieval  algorithm,  the 
KS77  relationship  is  usually  numerically  inverted  to  determine 
S  as  a  function  of  Tb  and  T s. 

For  a  given  emissivity,  Tb  linearly  varies  in  proportion  to 
Ts.  However,  Ts  (with  S )  influences  the  emissivity  itself, 
which  (for  both  V-  and  H-Pol  signals)  is  reduced  at  higher 
Ts.  Consequently,  for  salinities  higher  than  about  15  psu, 
Tb  first  increases  then  decreases  as  a  function  of  Ts ,  which 
makes  Ts  multivalued  with  respect  to  Tb  [a  line  parallel  to  the 
Ts-axis  intersects  an  isohaline  at  two  separate  values  of  Ts , 
Fig.  1(a)].  The  dashed  red  curve  identifies  the  maximum  Tb  for 
a  given  Ts.  In  contrast,  S  varies  uniquely  and  approximately 
inversely  (but  nonlinearly)  as  a  function  of  Tb ,  given  Ts. 
Hence,  given  Ts  and  an  observed  Tb ,  S  can  be  uniquely 
determined. 

The  SES  (plotted  in  Fig.  1(b)  on  axes  forming  an  oceanogra¬ 
pher’s  T—S  diagram)  uniquely  determines  density  D  as  a  non¬ 
linear  empirical  function  of  Ts  and  S.  Consequently,  if  Ts  and 
S  are  known,  D  can  be  determined.  Since  S  can  be  determined 
from  Tb  and  Ts  through  KS77,  the  two  relationships  can  be 
chained  together  to  retrieve  Z7,  given  Tb  and  Ts.  Hence,  D  may 
be  plotted  along  with  S  in  the  Tb-Ts  space.  Here,  density, 
which,  in  the  oceanographic  literature,  is  often  expressed 
in  units  of  Sigma-t  (crt,  the  density  anomaly  relative  to 


1000  kg  •  m-3,  e.g.,  10  at  corresponds  to  1010  kg  •  m*3)  is 
computed  as  a  function  of  temperature  and  pressure,  in  this 
case  at  the  surface. 

The  density  shown  in  relation  to  Tb  and  Ts  [Fig.  1(a)  (blue)] 
varies  like  5,  in  approximately  inverse  proportion  to  Tb.  It 
has  a  similar  but  more  pronounced  pattern  than  S  (red)  with 
stronger  negative  slopes,  i.e.,  Ts  is  multivalued  with  respect 
to  D  (a  line  parallel  to  the  Ts-axis  intersects  an  isopycnal  at 
two  separate  values  of  Ts)  and  reaches  a  maximum  (dashed 
blue  line)  at  various  Ts  values  throughout  its  range  (1000- 
1030  kg  •  m~3).  For  retrieving  S  and  D,  near  their  respective 
turning  points  (local  maxima),  the  algorithms  are  relatively 
insensitive  to  errors  in  Ts.  This  region  of  low  sensitivity 
appears  at  a  progressively  higher  Ts  for  density  and  salinity 
but  at  a  lower  Ts  for  D  than  for  S  (the  red  and  blue  dashed 
curves  are  offset,  with  respect  to  T s). 

The  error  sensitivity  of  S  and  D  is  explored  in  more  detail  in 
the  next  section,  but  we  can  immediately  conclude  that,  in  the 
presence  of  Ts  errors,  D  exhibits  low  sensitivity  to  Ts  at  lower 
temperatures  than  5  and  that,  for  both  S  and  D,  the  region  of 
low  sensitivity  appears  at  progressively  lower  temperatures  as 
they  increase  in  value.  The  sensitivity  to  Ts  errors  is  thus  lower 
for  S  in  relatively  fresh  ( S  ~  25  psu)  tropical  (Ts  ~  30  °C) 
waters  and  in  relatively  saline  ( S  ~  35  psu)  cool  temperate 
(Ts  15  °C)  waters.  For  density,  the  corresponding  areas  of 
low  Ts  sensitivity  occur  either  in  fresher  water  (S  ~  15  psu) 
or  at  lower  temperatures  (Ts  ~  23  °C)  with  respect  to  tropical 
waters,  whereas  the  differences  are  less  pronounced  in  cooler 
waters  (where  the  S  and  D  maxima  tend  to  converge). 

Given  the  unique  relationship  between  Tb  and  S  or  D,  we 
can  regard  Tb  (for  a  specified  incidence  angle  0)  as  being 
uniquely  determined  by  the  equation  of  state  [Fig.  1(b)]. 
However,  unlike  D,  which  monotonically  increases  as  S 
increases  and  Ts  decreases,  Tb  decreases  in  this  direction,  and 
Ts  is  multivalued  with  respect  to  T6,  on  the  more  saline  side  of 
the  diagram  ( S  >  20  psu);  a  line  parallel  to  the  vertical  Ts- axis 
intersects  a  (red)  Tb  isoline  at  two  separate  values  of  Ts. 

As  previously  mentioned  for  salinity,  the  relationships 
involving  Tb  and  D  were  computed  for  V-Pol  signals  at  an 
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(C)  (d) 


Fig.  2.  SSS  algorithm  sensitivities,  given  S  values  in  the  range  [1,5:5:40 
psuj  for  an  incidence  angle  of  22°,  vertical  polarization,  and  a  frequency  of 
1.4  GHz.  (a)  5  as  a  function  of  Ts  and  Tb.  Sensitivity  along  the  isohalines  of 
(b)  TbioTs  versus  T s ,  (c)  S  to  Tb  versus  Ts,  and  (d)  S  to  Ts  versus  Ts. 

incidence  angle  of  22°  relative  to  nadir.  The  plots  quantitatively 
differ  at  other  incidence  angles  but  show  the  same  qualitative 
features.  For  economy,  and  considering  that  STARRS  has  a 
single  V-Fol  mode,  we  exclusively  use  V-Pol  in  all  subsequent 
calculations  in  this  paper,  but  analogous  results  can  be  obtained 
for  other  polarization  modes. 

The  application  of  the  SSS  retrieval  algorithm,  which 
utilizes  the  KS77  relation  [replotted  in  Fig.  2(a)],  proceeds  by 
assuming  a  starting  (first  guess)  salinity  and  using  the  observed 
T s  value  to  compute  a  corresponding  Tb  for  a  given  instrument 
beam  incidence  angle  and  polarization.  Typically,  Ts  will 
either  be  measured  in  situ  or  remotely  from  satellite  or  airborne 
IR  or  microwave  radiometers,  with  appropriate  corrections 
for  atmospheric  effects,  to  ensure  that  it  is  accurate  at  the  sea 
surface.  The  computed  Tb  value  then  represents  the  brightness 
temperature  that  would  be  observed  by  the  instrument,  at  the  sea 
surface,  given  its  associated  incidence  angle  and  polarization 
state.  This  value  is  propagated  through  the  atmosphere  to  the 
remote  sensing  platform,  allowing  for  environmental  influences 
such  as  reflected  extraterrestrial  and  atmospheric  (O2)  radia¬ 
tion,  enhanced  emission  and  modified  reflection  due  to  a  rough¬ 
ened  sea  surface,  and  atmospheric  absorption.  The  resulting 
brightness  temperature  is  compared  with  the  observed  antenna 
temperature.  If  these  do  not  match  within  a  specified  tolerance, 
the  starting  salinity  is  adjusted  (based  on  a  gradient  descent 
algorithm),  and  the  process  iterates  until  convergence  is 
attained. 

The  development  of  an  alternative  algorithm  for  density  re¬ 
trieval  that  involves  a  direct  empirical  relationship  between  Ts, 
D,  and  Tb  for  a  given  S  and  0  could  also  be  considered.  Such 
an  algorithm  could  be  expressed  in  the  form  of  the  D  curves  in 
Fig.  1(a).  However,  we  are  aware  of  no  practical  technique  to 
directly  observe  density  in  the  ocean;  other  indirect  techniques 
such  as  measurement  of  refractive  index  are  less  accurate  and 
less  practical  than  the  conventional  methods  of  determining 
it  from  field  observations  of  Ts  and  S  using  a  well-calibrated 


conductivity  temperature  depth  (CTD)  profiler.  Direct 
determinations  of  seawater  density  in  a  laboratory,  e.g.,  through 
diluting  or  evaporating  seawater  and  measuring  corresponding 
masses  and  volumes,  are  also  difficult  to  accurately  perform 
in  comparison  with  laboratory  salinometer  measurement  of 
seawater  samples.  Given  that  the  practical  salinity  scale  [21] 
defines  S  in  terms  of  the  conductivity  of  a  standard  KC1  solution 
and  that  the  microwave  radiometers  respond  to  brightness 
temperatures  that  are  largely  determined  by  the  conductivity  of 
seawater,  the  approach  based  on  the  sequential  application  of 
the  empirical  KS77  and  SES  appears  to  be  the  most  practical. 

B.  Algorithm  Uncertainties 

For  the  chained  KS77  and  SES  density  retrieval  algorithm, 
uncertainties  will  arise  from  observational  biases  and  random 
errors  in  the  observed  Ts  and  T6,  which  are  input  to  the 
algorithm,  and  from  the  combined  effect  of  uncertainties  in 
the  empirical  relationships  for  determining  both  S  and  D. 
The  latter  uncertainties  could  arise  from  inaccuracies  in  the 
laboratory  data  as  well  as  from  errors  in  the  regression  model 
fits  used  to  determine  these  empirical  relationships. 

1)  Accuracy:  The  claimed  accuracy  of  KS77  based  on  var¬ 
ious  concentrations  of  standard  mean  ocean  water  is  0.3  K 
or  better  in  brightness  temperature.  While  the  KS77  relation 
was  determined  from  conductivity  measurements  using  a  mi¬ 
crowave  resonant  cavity  containing  samples  of  seawater  from 
Chesapeake  Bay  and  the  Pacific  Ocean  as  well  as  NaCl  solu¬ 
tions,  the  algorithm  has  not  yet  been  related  to  the  new  KC1 
conductivity  standard  upon  which  the  Practical  Salinity  Scale 
1978  [21]  is  based.  An  alternative  seawater  emissivity  model 
by  Ellison  et  al.  [8]  is  also  available.  This  employed  samples 
of  natural  and  synthetic  seawater  and  also  NaCl  solutions.  The 
latter  solutions  were  found  to  have  a  significantly  different 
permittivity  than  the  seawater  samples.  This  suggests  that  there 
might  be  difficulties  in  transferring  the  permittivity  algorithms 
to  the  new  KC1  standard.  The  Ellison  et  al.  emissivity  model 
reportedly  performs  better  than  KS77  at  high  microwave  fre¬ 
quencies  but  is  less  accurate  than  KS77  at  lower  microwave 
frequencies  (including  the  L-band).  Recent  laboratory  and  field 
experiments  by  Blanch  and  Aguasca  [22]  and  Wilson  etal.  [23] 
support  an  accuracy  of  ±0.1  K  for  KS77,  which  is  presently 
considered  the  de  facto  “standard”  for  work  at  L-,  S-,  and 
C-band.  These  are  the  most  useful  bands  for  surface  salinity, 
temperature,  and  roughness  retrievals  and  are  thus  most  relevant 
to  passive  remote  sensing  of  SSS. 

The  standard  error  associated  with  fitting  data  to  the  regres¬ 
sion  model  to  form  the  SES  is  given  as  3.6  x  10“3  kg  m-3 
in  [24].  The  accuracy  of  the  fitted  data  itself  was  not  reported, 
but  since  the  SES  is  the  virtual  standard  for  computing  sea¬ 
water  density  from  T  and  5,  we  assume  that  this  error  can 
be  neglected.  We  focus  instead  on  how  instrument  bias  and, 
particularly,  random  errors  for  Tb  and  Ts  propagate  through 
the  resulting  (chained)  nonlinear  retrieval  algorithm.  Given  the 
expected  observational  errors  in  remotely  sensing  SST  and  SSS, 
the  resulting  S  and  D  errors  turn  out  to  be  significantly  larger 
than  the  expected  errors  due  to  laboratory  data  inaccuracies  or 
regression  model  misfit  in  the  empirical  formulation  of  KS77 
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Fig.  3.  SSD  algorithm  sensitivities,  given  D  values  in  the  range 
[1000 :  5 :  1030  kg  m-3]  for  an  incidence  angle  of  22°,  vertical  polarization, 
and  a  frequency  of  1 .4  GHz.  (a)  D  as  a  function  of  Ts  and  Tb.  Sensitivity  along 
the  isopycnals  of  (b)  Tb  to  Ts  versus  Ts,  (c)  D  to  Tb  versus  Ts,  and  (d)  5  to 
Ts  versus  T s. 

and  the  SES.  Thus,  KS77  and  SES  model  errors  are  negligible 
compared  with  the  effects  of  present-day  measurement  error 
levels  in  Tb  and  Ts  when  deriving  errors  in  S  and  D. 

2)  Sensitivity:  Since  the  predicted  brightness  temperature 
Tb  is  (at  first  order)  proportional  to  the  SST  Ts  with  a  “con¬ 
stant”  of  proportionality  (the  emissivity),  being  itself  a  function 
of  temperature  and  salinity,  Tb  is  sensitively  and  nonlinearly 
related  to  Ts.  This  sensitivity  is  demonstrated  [Fig.  2(b)]  by 
the  absolute  value  of  the  derivative  of  Tb  with  respect  to  Ts 
versus  Ts  for  the  given  values  of  S ,  as  determined  by  KS77 
[Fig.  2(a)].  While  Figs.  1(a)  and  2(a)  present  Tb  as  a  function 
of  Ts,  Fig.  2(b)  presents  the  absolute  value  of  the  slope  (first 
derivative)  of  Tb  with  respect  to  Ts  versus  Ts.  Points  of 
zero  slope,  i.e.,  minimum  sensitivity,  correspond  to  the  turning 
points  traced  by  the  dotted  lines  in  Fig.  1(a).  These  points 
where  slope  changes  sign  from  positive  to  negative  at  higher 
Ts  values  appear  in  Fig.  2(b)  as  “reflections”  of  the  absolute 
value  of  the  slope  along  the  Ts-axis.  Points  of  maximum 
slope,  i.e.,  maximum  sensitivity,  which  appear  in  Fig.  1(a) 
where  the  parametric  curves  are  steepest,  correspond  to  the 
turning  points  (maxima)  of  the  curves  in  Fig.  2(b).  Fig.  2(b) 
thus  shows  that  the  minimization  and  maximization  of  the  Tb 
versus  Ts  slope  occur  within  certain  portions  of  the  Ts  range 
and  at  progressively  higher  temperatures  for  successively  lower 
salinities.  It  is  at  the  maxima  that  Ts  errors  will  most  adversely 
impact  the  determination  of  Tb. 

Corresponding  sensitivities  of  S  to  Tb  [Fig.  2(c)]  and  Ts 
[Fig.  2(d)]  as  a  function  of  Ts  exhibit  a  monotonically  changing 
negative  slope  as  a  function  of  increasing  Ts,  clearly  showing 
that  the  largest  SSS  retrieval  errors  due  to  instrumental  errors 
in  Tb  and  Ts  will  occur  at  the  colder  temperatures  and  the 
lower  salinities.  The  sensitivity  of  5  to  Ts  is  almost  zero 
at  the  highest  salinities  [Fig.  2(d)],  but  its  sensitivity  to  Tb 
remains  significant,  even  at  the  highest  Ts  and  S  values  (where 
it  asymptotes  to  about  1.6  psu/K).  From  the  perspective  of 
salinity  retrieval,  we  note  that  high  sensitivity  of  5  to  Ts  is 
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Fig.  4.  Idealized  horizontal  transects  of  (a)  SST  and  (b)  Tb  with  given  rms 
error  levels  and  the  corresponding  (c)  SSD  and  (d)  SSS  with  computed  95% 
error  bounds.  The  Ts  and  Tb  profiles  were  arbitrarily  chosen,  so  that  horizontal 
variations  in  Ts  and  S  mutually  reinforce  or  compensate  that  of  D  at  different 
locations.  An  incidence  angle  of  22°  was  assumed. 

a  disadvantage  to  the  extent  that  errors  in  Ts  will  adversely 
impact  the  retrieval,  whereas  low  or  even  zero  sensitivity  is 
advantageous.  Similarly,  high  sensitivity  of  S  to  Tb  is  a  disad¬ 
vantage,  because  it  amplifies  the  effects  of  Tb  errors,  whereas 
low  sensitivity  is  advantageous.  In  this  case,  zero  sensitivity 
would  be  undesirable,  because  a  complete  lack  of  dependence 
of  S  on  Tb  would  make  the  retrieval  of  S  from  Tb  impossible! 

For  the  derivation  of  density  from  the  combined  empirical 
relations  for  SSS  and  D  [Fig.  3(a)],  the  patterns  of  sensitivity 
to  the  observables  Tb  and  Ts  are  similar  to  those  of  5,  as 
previously  discussed.  The  sensitivity  of  Tb  to  Ts,  given  D 
[Fig.  3(b)],  although  similar  to  that  given  S  [Fig.  2(b)],  has 
both  its  turning  and  reflection  points  shifted  toward  colder  Ts 
values  and  a  maximum  occurring  for  the  least  dense  water 
(1000  kg  •  m~3)  at  the  high-temperature  end.  This  is  due  to 
the  slopes  of  D  tending  to  be  more  negative  than  those  of  S 
at  higher  Ts,  as  revealed  in  Fig.  1(a)  (where  the  turning  points 
for  D  traced  by  the  dotted  lines  fall  to  the  left  of  those  of  S). 
The  sensitivities  of  D  to  Tb  and  Ts  [Fig.  3(c)  and  (d)],  as  a 
function  of  Ts,  all  show  similar  patterns  to  those  of  S. 

To  illustrate  how  random  observational  errors  in  Ts  and  Tb 
impact  the  retrieval  of  S  and  D  under  a  representative  range 
of  T  and  Tb  values,  we  define  an  idealized  spatial  transect  of 
Ts  and  Tb ,  with  an  arbitrary  distance  scale  and  origin  ( d  =  0) 
at  the  center  of  the  transect,  that  could  represent  a  strong 
continental  shelf  scale  frontal  contrast  (Fig.  4).  This  might 
be  encountered  in  wintertime  shelf  circulation  when  a  warm 
western  boundary  current,  such  as  the  Gulf  Stream,  impinges  on 
the  outer  continental  shelf.  In  this  case,  Ts  [Fig.  4(a)]  is  chosen 
to  increase  monotonically  seaward,  whereas  Tb  similarly,  but 
arbitrarily,  increases  in  the  same  direction  and  finally  decreases 
at  the  seaward  end.  These  “observed”  Ts  and  Tb  values  were 
used  to  perform  a  simulated  salinity  and  density  retrieval  in  the 
presence  of  representative  random  (root  mean  square,  rms)  er¬ 
rors  in  the  raw  Ts  and  Tb  observations  of  0.25  °C  and  0.5  K, 
respectively.  Although  arbitrarily  chosen,  these  values  approx¬ 
imate  those  estimated  for  30-day  averaged  SMOS  data  [see 
Section  II-C  and  Table  1(b)].  The  resulting  mean  D  distribution 
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TABLE  I 

Input  Sampling  Parameters,  and  SST,  WS, 
and  Tb  Errors  for  SSD  Retrieval 


(a) 


Instrument 

IK  SST 

Stdev 

K 

L-band 

Tb 

Stdev 

K 

l-band 
Integration 
Time  s 

Averaging 

Period 

SST  (& 
WS) 

A* 

Tb 

A* 

STARRS 

0.70 

0.21 

24.0 

1  min 

1 

2.5 

SMOS 

0.70 

4.00 

1.2 

30  day 

10 

50 

Aquarius 

0.70 

0.06 

12.0 

28  day 

4 

4 

(b) 

Nominal 

Ts  noise 

Tb  Noise 

Wind  Spd 

Wind 

Total  Tb 

Piiel  Size 

Stdev 

Stdev 

Stdev 

Noise 

Nobe 

km 

K 

K 

m/s 

Stdev 

Stdev 

K 

K 

STARRS 

5 

0.70 

0.13 

1.5 

0.21 

0.25 

SMOS 

50 

0.22 

0.57 

1.5 

0.03 

0.57 

Aquarius 

80 

•J5 

0.03 

1.5 

0.06 

0.07 

(C) 

Incidence 

Ts  bias 

Tb  Bias 

Wind  Spd 

Wind 

Total  Tb 

Angle 

Rms 

Rms 

Bias  Rms 

Tb  Bias 

Bias  Rms 

deg 

K 

K 

m/s 

K 

K 

STARRS 

22 

0.10 

0.71 

0.5 

0.07 

0.71 

SMOS 

40 

0.20 

0.40 

0.5 

0.03 

0.40 

Aquarius 

34 

0.20 

0.10 

0.5 

0.04 

0.11 

[Fig.  4(c)]  shows  a  substantial  decrease  offshore  that  is  driven 
by  the  Ts  increase  (with  warmer  water  of  the  same  salinity 
being  less  dense),  but  it  starts  to  increase  at  the  outer  end 
while  Ts  merely  levels  off.  Reference  to  the  retrieved  salinity 
[Fig.  4(d)]  shows  that  the  density  trend  reversal  was  associated 
with  an  increase  in  S  at  the  outer  shelf  end  of  the  transect  (near 
d  =  40  km).  Confidence  bands  (at  the  95%  level)  given  for 
rms  error  levels  corresponding  to  raw  and  (three-  and  five-point 
box  car)  filtered  data  are  shown  for  the  input  T  and  Tb  values 
and  the  output  D  and  S  values.  While  these  bands  are  propor¬ 
tionally  small  and  uniform  across  shelf  for  Ts  and  Tb ,  they 
are  proportionately  wider  but  decrease  in  width  across  shelf 
for  D  and  S.  The  reduction  in  rms  error  in  the  outer  shelf 
is  due  to  the  decreased  sensitivity  of  D  and  S  to  Ts  and  Tb 
errors  experienced  in  warmer  waters  [Figs.  2(c)  and  (d)  and 
3(c)  and  (d)]. 

Both  the  S  and  D  retrieval  errors  were  computed  using  two 
different  methods.  The  first  assumed  that  both  the  input  and 
output  errors  were  normally  distributed  and  defined  the  95% 
confidence  limits  as  being  twice  the  standard  deviation  of  the 
series.  The  retrievals  were  then  performed  on  both  the  mean 
values  and  the  upper  and  lower  bounds  of  Ts  and  Tb.  In  this 
way,  the  error  limits  were  propagated  through  the  retrieval 
algorithm  by  applying  the  same  equations  to  the  limits  that 
were  used  to  transform  the  mean  values.  The  second  method 
employed  Monte  Carlo  simulation  to  deduce  the  rms  errors. 
This  was  done  by  generating  a  large  number  of  realizations  of 
the  Ts  and  Tb  space  series  using  error  time  series  produced 
by  a  random  normal  number  generator.  The  retrievals  were 
then  performed  on  each  realization  of  the  ensemble,  and  the 


results  were  statistically  analyzed  to  determine  the  mean  and 
rms  errors  of  D  and  S.  Since  both  methods  produced  virtually 
identical  results,  the  first  method  was  subsequently  employed, 
because  it  was  computationally  faster. 

C.  Density  Retrieval  Errors 

Representative  SSD  retrieval  errors,  both  systematic  and 
random,  are  now  estimated  for  selected  sampling  configurations 
and  parameters  of  SMOS  and  Aquarius.  These  are  needed  to 
evaluate  potential  applications  of  the  satellite  sensors  and  to 
determine  the  utility  of  STARRS  for  testing  satellite  retrieval 
algorithms.  The  diverse  technologies,  sampling  capabilities, 
and  operational  configurations  of  these  instruments  do  not 
allow  a  comprehensive  analysis.  This  would  require  the  appli¬ 
cation  of  various  simulation  methods,  e.g.,  [25]— [29],  and  is 
outside  our  scope.  Instead,  we  select  representative  sampling 
parameters  (Table  I)  for  typical  operating  conditions  of  each 
system.  The  intent  is  not  to  compare  and  contrast  sensor  per¬ 
formance,  which  varies  over  a  wide  parameter  range  and  can 
be  optimized  for  specific  purposes,  but  to  indicate  expected 
error  levels  when  retrieving  SSD  under  conditions  that  are 
typical  for  each  instrument.  Although  several  environmental 
corrections  are  needed  to  minimize  noise  and  bias  in  SSS  and 
SSD  retrievals  (see  Section  II-A),  we  focus  on  SST  and  wind- 
induced  sea  surface  roughness  (SSR),  which  exert  the  most 
influence  on  the  retrieval  errors.  STARRS  Ts  and  Tb  errors 
were  estimated  from  laboratory  and  field  studies  ([30],  whereas 
representative  Ts ,  Tb ,  and  associated  wind-induced  noise  and 
bias  values  for  SMOS  and  Aquarius  were  obtained  from  the 
literature  (e.g.,  [26],  [27],  [31],  and  [32]).  We  discuss  random 
errors  first  and  then  show  how  they  combine  with  systematic 
errors  to  determine  the  total  SSD  error  due  to  Ts  and  Tb. 

Estimated  single-sample  errors,  which  are  expressed  as  stan¬ 
dard  deviations  Sdv  for  Ts  and  Tb  Table  1(a)  (columns  1  and 
2)  were  determined  from  laboratory  and  field  measurements. 
The  SST  values  are  the  Sdv' s  of  the  residual  errors  of  obser¬ 
vations  from  the  onboard  IR  radiometer  (for  STARRS)  and 
from  the  National  Oceanic  and  Atmospheric  Administration 
(NOAA)  Advanced  Very  High  Resolution  Radiometer  data 
(for  the  satellites),  compared  with  in  situ  observations,  after 
atmospheric  water  vapor  correction.  Tb  values  correspond  to 
the  effective  noise  equivalent  delta  temperature  (or  radiometric 
sensitivity)  of  each  L-band  radiometer  for  the  specified  integra¬ 
tion  time  (column  4).  For  the  random  error  analysis  [Table  1(b)], 
the  single-sample  data  acquired  at  the  nominal  pixel  sizes 
(column  2,  given  for  illustration  but  not  used  in  the  calculations) 
are  assumed  to  be  temporally  averaged  over  a  given  period 
[Table  1(a)  (column  5)].  For  STARRS,  this  implies  spatial  aver¬ 
aging  over  a  5-km  distance  along  track  in  a  single  overpass  (the 
actual  sampling  times  and  pixel  sizes  are  ~2  s  and  1  km,  respec¬ 
tively;  see  Appendix  A-l).  We  could  average  successive  daily 
overflights  of  STARRS  over,  for  example,  a  ten-day  period, 
but  salinity  distributions  in  coastal  areas,  where  it  is  typically 
flown,  vary  significantly  on  daily  time  scales,  so  aliasing  errors 
would  be  large.  For  satellite  observations  over  the  deep  ocean, 
which  vary  on  much  longer  time  scales,  averaging  of  data  from 
particular  pixels  during  successive  repeat  orbits  is  assumed.  For 
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SMOS,  which  is  an  imaging  interferometer,  successive  multi¬ 
angle  views  acquired  along  the  “dwell  line”  provide  additional 
samples  for  temporal  averaging  [27].  Hence,  the  corresponding 
number  of  degrees  of  freedom  JVdf  (or  the  effective  number 
of  independent  samples)  for  Tb  [Table  1(a)  (column  6)]  is  the 
least  for  STARRS  and  the  most  for  SMOS.  SST  and  wind 
observations  are  assumed  to  be  acquired  near  each  L-band  pixel 
only  once  per  orbit,  so  their  ATdf  [Table  1(a)  (column  4)]  is 
generally  smaller  than  that  for  Tb.  For  SMOS  and  Aquarius, 
temporal  averaging  is  more  efficient  than  spatial  averaging  due 
to  spatially  correlated  environmental  correction  and  instrument 
errors  [27].  While  the  actual  pixel  sizes  and  incidence  angles  for 
each  instrument  vary  over  a  significant  range  (e.g.,  32-100  km 
and  0°  to  65°  for  SMOS  [5]),  only  nominal  midrange  values  are 
used  here.  Although  the  data  from  both  the  SMOS  and  Aquarius 
missions  will  be  spatially  averaged  at  ~300-km  resolution  to 
produce  monthly  global  salinity  products,  with  a  target  pre¬ 
cision  in  the  range  of  0.1  -0.2  psu,  we  assume,  for  simplicity 
and  to  preserve  resolution,  that  no  spatial  averaging  is  per¬ 
formed  (with  the  exception  of  STARRS,  as  previously  noted). 
To  account  for  noise  reduction  by  temporal  averaging,  we 
assume  independent  Gaussian  random  samples  and  reduce  the 
sampling  Sdv  [Table  1(a)]  by  a  multiplicative  factor  of 
to  determine  effective  noise  levels  Table  1(b)  (columns  3 
and  4)  for  IR  Ts  and  L-band  Tb ,  respectively.  For  the  L-band, 
there  is  an  added  contribution  from  (primarily  wind-induced) 
SSR.  The  wind  speed  random  error  of  1 .5  m  •  s-1  (column  5)  is 
representative  of  available  buoy  data  and  satellite-observed  or 
model-analyzed  ocean  wind  products.  Using  this  value  with  the 
incidence  angles  [Table  1(c)  (column  1)]  in  the  Camps  et  al. 
[33]  empirical  formula  (for  V-Pol)  yields  the  wind  impact 
on  rms  Tb  noise  levels  [Table  1(b)  (column  6)].  These  are 
added  to  the  instrument  Tb  noise  (using  root  sum  of  squares 
addition)  to  yield  the  total  Tb  random  noise  level  [Table  1(b) 
(column  7)].  The  Ts  (column  3)  and  Tb  (column  7)  noise 
levels  are  separately  input  into  the  algorithm  to  determine  their 
impact  on  random  errors  in  SSS  and  SSD  retrievals  through 
the  combined  effect  of  KS77  and  the  SES.  For  each  T -S  pair 
in  a  preassigned  matrix  spanning  their  observational  range,  the 
sensitivity  of  SSS  and  SSD  (Figs.  2  and  3)  to  errors  in  Ts  and 
Tb  is  multiplied  by  the  corresponding  Ts  and  Tb  noise  levels 
Table  1(b)  (columns  3  and  7)  to  produce  a  map  of  SSD  random 
errors  in  the  T -S  space  [Fig.  5(a),  (c),  and  (e)]. 

In  contrast  to  random  errors,  which  are  reduced  by  averaging, 
bias  estimates  are  not  reduced.  Being  thus  simpler  to  statisti¬ 
cally  treat,  they  are  harder  to  estimate  in  practice,  because  of  the 
difficulty  of  establishing  absolute  reference  values  in  the  field. 
The  biases  for  Ts  [Table  1(c)  (column  3)]  were  estimated  from 
published  comparisons  of  water- vapor-corrected  SST  product 
values  with  in  situ  buoy  data,  e.g.,  [34]  and  [35],  whereas  the 
instrument  Tb  biases  (column  4)  were  determined  from  actual 
laboratory  data  (STARRS  [30])  and  from  simulation  or  pro¬ 
totype  tests  (SMOS  and  Aquarius  [31],  [32])  giving  expected 
calibration  drift  over  one  month.  The  wind  bias  (column  5)  was 
estimated  from  comparisons  of  various  satellite  wind  sensor 
products  and  buoy  winds,  e.g.,  [36].  The  instrument  and  wind- 
induced  Tb  rms  biases  and  noise  levels  were  finally  added 
together  (in  the  rms  sense),  to  determine  the  total  Tb  error 


(column  7)  for  input  into  the  retrieval  algorithm,  along  with  the 
corresponding  Ts  errors  (column  3). 

The  resulting  SSD  error  maps  (Fig.  5)  exhibit  a  broad  ellipti¬ 
cal  minimum  stretching  between  the  high  T  and  5-axis  values, 
with  a  focus  mostly  near  high  T.  The  errors  generally  increase 
as  SST  and  SSS  values  decrease  and  ultimately  reach  levels 
that  are  unacceptably  large  for  any  conceivable  application 
(in  excess  of  3  kg  •  m“3).  The  total  (random  plus  systematic) 
error  shown  in  the  right-hand  panels  is  minimized  for  warm 
SSTs  and  midrange  SSS,  and  is  low  (~0.25)  for  Aquarius, 
due  mainly  to  its  exceptional  Tb  stability  and  noise  perfor¬ 
mance.  It  is  marginal  (exceeding  1  kg  m-3)  in  the  midranges 
for  STARRS  and  SMOS.  The  random  errors  isolated  in  the 
left-hand  panels  follow  a  similar  distribution  but  are  signifi¬ 
cantly  lower  than  the  total  errors  for  SMOS,  and  particularly 
STARRS,  which  has  a  relatively  good  noise  performance  but 
poorer  calibration  stability.  The  SSD  error  variation  (black 
contours)  is  approximately  orthogonal  to  that  of  density  (white 
contours),  so  the  two  quantities  are  virtually  uncorrelated.  In 
any  case,  the  density  errors  are  on  the  order  of  1  :  1 000  kg  m  3 
or  0.1%.  However,  the  errors  expressed  in  units  of  ot  as  a 
percentage  of  the  range  of  density  variation  of  the  world’s  ocean 
~35  <7t(1035  kg  •  m“3  -  1000)  are  more  pessimistic,  being 
~1%  at  best  and  rapidly  climbing  above  10%  at  SSS  levels  of 
below  about  15  psu,  which  is  typical  of  estuarine  plumes.  The 
SSD  error  levels  for  all  three  systems  ~0.5  kg  m~3  (or  0.5  at) 
for  the  T  and  S  values  typical  of  temperate  and  tropical  climates 
appear  adequate  for  determining  density  gradients,  at  least  in 
more  saline  areas,  and  where  density  fields  vary  significantly 
on  scales  exceeding  100  km. 

III.  Density  Variation 

Based  on  the  preceding  idealized  SSD  error  analysis  and 
the  more  realistic  cases  of  STARRS,  SMOS,  and  Aquarius, 
we  now  consider  how  SST  and  SSS  variations  and  errors  may 
compensate  or  reinforce  one  another  to  influence  retrieved  SSD 
and  its  errors.  We  then  evaluate  the  potential  to  determine  SSD 
gradients  with  sufficient  precision  to  determine  near-surface 
geostrophic  shear  from  SST  and  SSS  mapped  by  the  IR  and 
microwave  satellites. 

A.  Compensation  and  Reinforcement 

The  somewhat  arbitrary  distribution  of  Ts  and  Tb  [Fig.  4(a) 
and  (b)]  resulted  in  a  retrieved  S  [Fig.  4(d)]  that  was  largely 
anticorrelated  with  Ts  (except  in  the  outer  shelf  where  Ts 
leveled  off).  This  represents  an  interesting  case  of  mutual 
reinforcement  in  which  the  T  and  S  variations  both  contribute 
to  a  decrease  in  D  (as  T  increases  and  S  decreases).  In  this 
situation,  D  variations  will  always  be  large  when  (both)  Ts  and 
S  variations  are  large.  Consequently,  frontal  systems  may  be 
best  detected  by  identifying  regions  of  strong  density  variation. 
The  reinforcement  case  may  occur  in  summer  when  warm 
water  near  shore  contributes  to  low  near-shore  densities  that  are 
reinforced  by  runoff  from  large  rivers  with  persistent  annual 
flow  (or  little  seasonal  variation).  This  could  be  illustrated  by 
the  case  shown  in  Fig.  4  if  the  coast  is  considered  to  lie  on 
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Fig.  5.  (a),  (c)  and  (e)  Random  and  (b),  (d)  and  (0  total  (i.e.,  random  plus  systematic)  SSD  error  (grayscale  with  black  contour  lines)  as  a  function  of  SST  and 

SSS  for  representative  configurations  of  (a)  and  (b)  STARRS,  (c)  and  (d)  SMOS,  and  (e)  and  (f)  Aquarius.  The  white  contour  lines  with  markers  represent  SSD 
ranging  from  1000  to  1030  (left  to  right)  in  5-kg  m-3  steps  [as  in  Fig.  1(b)]. 


the  right-hand  side  (RHS),  rather  than  on  the  left,  as  in  the 
preceding  discussion. 

The  opposite  case  of  compensation  is  also  surprisingly  com¬ 
mon  in  oceanographic  settings.  For  example,  in  temperate 
regions  during  winter,  cooler  water  near  the  coast  may  coincide 
with  freshwater  associated  with  winter  runoff.  Cases  of  com¬ 
pensation  also  occur  in  frontal  systems  that  are  equilibrated  by 
strong  vertical  mixing  processes  operating  at  the  frontal  inter¬ 
face  (e.g.,  salt  fingering).  In  the  Middle  Atlantic  Bight  of  the 
U.S.  east  coast,  cool  freshwater  likely  derived  from  the  ice  melt 


in  the  Labrador  Sea  (forming  the  so-called  “Cold  PooP;  see 
[12]  for  sample  sections)  may  lie  shoreward  of  and/or  under¬ 
neath  the  warm  saline  slope  water  derived  in  part  from  the  Gulf 
Stream.  In  such  cases  (e.g.,  Fig.  6),  Ts  and  S  both  increase  sea¬ 
ward,  so  their  variations  are  correlated  and  they  have  a  mutually 
compensating  effect  on  the  resulting  D  variation.  Consequently, 
the  density  contrast  may  be  weak  or  even  absent,  whereas  Ts 
and  S  variations  are  strong,  and  fronts  may  not  be  detected  in 
density  retrievals,  even  though  they  are  strongly  evident  in  Ts 
and  S  retrievals.  In  Fig.  6,  we  dispensed  with  the  error  bounds. 
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Weak  Compensation  Case 


Fig.  6.  Idealized  horizontal  transects  of  (a)  SST  and  (b)  SSS,  and  the 
corresponding  (c)  SSD  and  (d)  Tb.  The  Ts  and  S  profiles  were  chosen,  so 
that  their  horizontal  variations  mutually  and  weakly  compensate  in  their  effect 
on  D.  The  corresponding  Tb  for  incidence  angles  of  7°,  22°,  and  37°  are  also 
shown  in  (d).  Error  bounds  for  D  in  this  case  are  shown  in  Fig.  9(c). 


in  favor  of  showing  how  Tb  varies  with  incidence  angle. 
However,  the  uncertainty  in  the  D  variation  is  illustrated  for  the 
arbitrary  and  weak  compensation  cases  in  Figs.  9  and  10.  For 
the  latter  case,  the  across-shelf  density  gradient  is  not  resolved. 


B.  Density  Compensation  Index 

Given  the  importance  of  density  reinforcement  and  compen¬ 
sation  processes  in  physical  oceanography  and  the  emerging 
possibility  of  determining  density  from  salinity  mapping  satel¬ 
lites,  it  is  instructive  to  employ  an  index  that  can  be  mapped 
in  a  survey  area  to  help  determine  the  potential  dynamical 
significance  of  ocean  fronts  associated  with  strong  temperature 
and  salinity  contrasts.  Apart  from  knowing  whether  T  and 
S  are  mutually  reinforcing  or  compensating  with  respect  to 
density,  it  is  also  of  interest  to  determine  which  of  T  and  S 
play  the  dominant  role  in  determining  D.  This  can  be  useful 
in  developing  numerical  hydrodynamic  models  to  simulate 
likely  contributions  to  model  density  fields  from  such  diverse 
factors  as  freshwater  runoff,  solar  heating,  and  ice  melt.  To 
serve  both  these  purposes,  we  employ  a  vector  index  (the 
“Turner  vector,”  presented  in  Appendix  B),  which  is  a  simple 
extension  of  the  Turner  angle  defined  by  Ruddick  [18]  for 
studying  double  diffusive  activity.  This  angle  is  defined  in  terms 
of  the  density  ratio  used  by  oceanographers  to  study  ocean 
stratification  and  fine  structure.  Johnson  [19]  adapted  the  Turner 
angle  to  study  density  compensation  in  horizontal  T  and  S 
gradients.  He  projected  the  local  horizontal  salinity  gradient 
onto  the  direction  of  the  local  horizontal  surface  temperature 
gradient  (which,  in  the  deep  ocean,  is  mostly  equatorward). 
This  provided  the  reference  direction  for  computing  the  “hor¬ 
izontal”  Turner  angle,  enabling  large  areas  of  the  world’s  ocean 
with  density  compensating  salinity  gradients  to  be  identified. 
In  our  coastal  work,  the  across-shelf  direction  is  a  convenient 
reference  direction,  which  closely  approximates  the  directions 
of  the  prevailing  T,  5,  and,  hence,  D  gradients  through  the 
high-aspect-ratio  Plata  plume. 

The  Turner  vector  can  be  used  to  indicate  either  density 
reinforcement  or  compensation  by  T  and  S  while  taking  into 


account  the  relative  dominance  and  absolute  significance  of  T 
and  S  variations  in  their  effect  on  D.  To  illustrate  its  application 
to  various  idealized  oceanographic  settings,  we  consider  a 
variety  of  cases  including  strong  reinforcement  with  salinity 
dominance  [Fig.  7(a)];  weak  compensation  with  S  dominance 
[Fig.  7(b)];  the  arbitrary  variation  switching  from  reinforce¬ 
ment  to  compensation  with  S  dominating,  except  during  the 
transition  when  T  dominates  for  a  short  distance  [Fig.  7(c)]; 
and  finally  the  case  of  a  linear  variation  in  density  across 
shelf,  which  results  in  reinforcement  mostly  dominated  by 
S  variations  [Fig.  7(d)].  The  cases  of  weak  compensation 
and  arbitrary  variation  are  those  presented  in  Figs.  4  and  6, 
respectively.  The  strong  reinforcement  case  was  created  from 
arbitrary  but  opposing  T  and  S  gradients.  The  linear  case  was 
created  from  an  arbitrary  nonlinear  across-shelf  T  and  a  linear 
D  pattern  with  the  corresponding  Tb  (and  S)  computed  using 
the  SES  and  KS77  algorithms  in  the  forward  (instead  of  the 
usual  inverse)  computation  modes.  This  somewhat  artificial 
reinforcement  case  was  created  to  illustrate  how  a  linear  D 
variation  in  space  could  (and  in  fact  must)  result  from  nonlinear 
T  and  S  variations  due  to  the  inherent  nonlinearities  of  the 
algorithm.  The  perfect  compensation  case  ( dD/dx )  =  0  is  a 
special  case  of  a  linear  D  variation,  with  zero  across-shelf 
gradient.  These  two  cases  also  illustrate  how  quite  different 
patterns  of  T  and  S  reinforcement  or  compensation  can  produce 
density  distributions  that  are  similar  in  character,  i.e.,  D  profiles 
are  not  uniquely  determined  by  T  and  S. 

The  Turner  vector  values  for  the  cases  shown  in  Fig.  7  may 
be  displayed  more  compactly  in  a  polar  diagram  (Fig.  8),  which 
is  a  variant  of  the  approach  used  by  Ruddick  [18],  which  is 
explained  in  Appendix  B.  This  method  of  presentation  displays 
variations  of  the  Turner  angle  and  magnitude  for  different 
transects  in  a  compact  diagram  but  does  not  explicitly  reveal 
the  pattern  of  across-shelf  variation  evident  in  Fig.  7. 

IV.  Horizontal  Gradient  and  Geostrophy 

In  the  previous  sections,  we  considered  how  SSS  and  SSD 
may  be  derived  using  existing  salinity  retrieval  algorithms, 
combined  with  the  SES.  We  also  examined  associated  retrieval 
uncertainties  and  considered  the  nonlinear  influence  of  various 
spatial  distributions  of  Ts  and  Tb  on  the  spatial  distribution 
of  D.  The  resulting  density  gradients  may  be  weak  or  strong, 
depending  on  whether  Ts  and  S  are  either  weak  or  mutu¬ 
ally  compensating,  or  else  mutually  reinforcing.  These  have 
significant  impact  on  the  dynamics  of  ocean  circulation.  One 
significant  impact,  which  is  widely  recognized  by  physical 
oceanographers,  is  its  effect  on  horizontal  pressure  gradients 
and,  hence,  on  buoyancy-driven  currents.  Geostrophy  theory 
depends  on  a  simplified  momentum  balance,  in  which  the 
Coriolis  and  horizontal  pressure  gradient  terms  cancel,  whereas 
the  other  terms  are  negligible.  Widely  used  in  classical  mete¬ 
orology  and  oceanography  to  determine  mean  ocean  currents 
and  winds,  it  has  recently  come  to  the  forefront  again,  with 
the  advent  of  high-performance  spaceborne  altimeters  (e.g., 
TOPEX  and  Jason).  These  can  resolve  SSH  with  an  accuracy 
of  about  5  cm  and  a  precision  of  2  cm  or  better.  They  can 
be  used  to  determine  sea  surface  slope  variations  with  respect 
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Fig.  7.  Turner  vector  plots  shown  in  Cartesian  component  form  as  a  function  of  across-shelf  location.  The  cases  shown  are  idealized  transects  exhibiting 
(a)  reinforcement,  (b)  weak  (under  )  compensation,  (c)  arbitrary  density  variation,  and  (d)  linear  density  variation.  Cases  (b)  and  (c)  correspond  to  the  profiles 
shown  in  Figs.  6  and  4,  respectively. 


Comp 

dp/dx=0  Turner  Angle 


Comp 

Fig.  8.  Turner  vector  plots  shown  in  polar  form,  corresponding  to  the  idealized 
transects  shown  in  Fig.  7,  i.e.,  (*)  reinforcement,  (+)  weak  (under-)  compensa¬ 
tion,  (o)  arbitrary  density  variation,  and  (x)  linear  density  variation. 

to  the  mean  ocean  surface  topography,  from  which  surface 
geostrophic  currents  can  be  determined.  Along  with  accurate 
geoid  observations  from  gravity  missions  such  as  Champ  and 
GRACE,  this  makes  possible  the  determination  of  absolute 
(i.e..  Earth  relative)  mean  and  fluctuating  ocean  currents  due 
to  horizontal  pressure  gradients  expressed  at  the  surface. 

The  possibility  of  mapping  SSD  from  space  raises  the 
prospect  of  computing  the  vertical  shear  of  the  geostrophic 
currents  at  the  surface  to  give  an  indication  of  subsurface 
current  profiles  driven  by  horizontal  pressure  gradients  arising 
from  surface  and  internal  density  differences.  If  we  assume 
that  horizontal  density  gradients  dD(x)/dx  are  independent 
of  depth  (at  least,  within  a  depth  range  of  interest),  while  the 


density  linearly  increases  with  depth,  we  obtain  a  formula  for 
the  total  geostrophic  current  Vg  at  depth  2  (see  Appendix  C  for 
details  of  its  derivation  and  notation),  i.e.. 


_  9_  d£  _  gz  dD(x) 
f  dx  fD(x)  dx 


(1) 


Here,  g  and  /  are  the  acceleration  due  to  gravity  and  the 
Coriolis  constant,  respectively,  and  f  is  the  sea  surface  elevation 
relative  to  the  mean  sea  level.  The  two  right-side  terms,  which 
we  denote  as  Ve  and  V i,  result  from  pressure  gradients  driven  by 
the  external  sea  surface  slope  and  internal  density  differences, 
respectively.  Ve  can  be  obtained  by  deriving  SSH  gradients 
from  satellite  altimeters,  whereas  Vi  can  be  obtained,  as  demon¬ 
strated  in  this  paper,  by  deriving  horizontal  density  gradients 
from  L-band  radiometers  using  a  variant  of  the  “thermal  wind" 
equation.  At  the  surface,  Vi  =  0,  so  the  surface  geostrophic 
velocity  V^(0)  =  Ve ,  whereas,  at  depth  2,  Vi(z )  represents  the 
deviation  of  the  geostrophic  velocity  from  that  at  the  surface. 
In  meteorology  and  oceanography,  the  term  “velocity  shear" 
is  commonly  used  to  describe  this  deviation,  but  we  prefer  to 
apply  this  term  to  the  vertical  velocity  gradient  dVg/dz. 

By  assuming  that  dD(x)/dx  is  constant  throughout  a  near¬ 
surface  layer  of  specified  depth,  we  can  project  the  vertical 
shear  at  the  surface  downward  to  determine  the  geostrophic 
current  throughout  the  layer.  This  assumption  may  be  expected 
to  hold  in  a  variety  of  situations  such  as  in  continental  shelf 
fronts  during  winter  when  the  large-scale  stratification  is  not 
perturbed  by  the  seasonal  pycnocline  and  along  the  margins 
of  the  deep  ocean  gyres  away  from  intense  western  boundary 
currents,  which  may  be  bounded  by  nonlinear  meanders  and 
mesoscale  eddies.  We  will  show  later  that  this  assumption  is 
justified  in  the  large-scale  wintertime  Plata  river  plume.  The 
assumption,  together  with  the  results  of  a  dimensional  analysis 
demonstrating  that  the  across-plume  momentum  balance  is 
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Fig.  9.  Derivation  of  geostrophic  current  at  the  base  of  the  surface  layer  with 
95%  error  bounds  for  simulated  raw  and  filtered  data  for  the  case  of  arbitrary 
horizontal  SST  and  SSS  variation.  The  raw  rms  errors  for  Ts  and  Tb  were 
0.25  °C  and  0.5  K,  respectively.  The  box  car  filters  were  one,  three,  and 
five  points  long.  The  transects  shown  are  for  surface  (a)  density,  (b)  horizontal 
density  gradient,  (c)  geostrophic  shear,  and  (d)  geostrophic  velocity  at  depth 
(zo).  relative  to  the  surface.  This  case  corresponds  to  the  arbitrary  SSS  and 
SST  transects  shown  in  Fig.  4. 

geostrophic,  can  he  used  to  estimate  the  subsurface  geostrophic 
current  on  the  continental  shelf. 

If,  for  the  moment,  we  assume  that  the  surface  geostrophic 
velocity  Ve  =  0,  then  the  geostrophic  current  at  the  base  of  a 
near-surface  layer  of  thickness  H  is  just  the  internal  velocity 
given  by  (1),  i.e., 

v*<">  -  - -jm '  ^  121 

This  calculation  may  be  illustrated  by  using  an  idealized 
case,  which  also  allows  the  resulting  error  bounds  on  the 
geostrophic  velocity  to  be  determined.  For  the  arbitrary  case 
presented  earlier  [Fig.  9(a)],  the  horizontal  density  gradient 
[Fig.  9(b)]  is  negative  near  the  center  of  the  transect,  with 
smaller  errors  in  both  D  and  dD/dx  due  to  the  warmer  temper¬ 
atures  at  the  outer  shelf  end  [Fig.  4(a)].  The  random  errors  are 
amplified  by  the  differencing  operation;  hence,  spatial  filtering 
is  necessary.  Since  bias  errors  due  to  instrument  drift  will  be 
largely,  if  not  completely,  eliminated  by  the  differencing  oper¬ 
ation,  we  have  chosen  not  to  simulate  them  here.  The  resulting 
geostrophic  shear  (or  vertical  velocity  gradient)  [Fig.  9(c)]  and 
geostrophic  current  (relative  to  the  surface)  at  the  base  of  a 
200-m-deep  layer  [Fig.  9(d)]  show  that,  for  the  specified  Ts 
and  Tb  error  levels,  and  modest  (three-  or  five-point  box  car) 
filtering,  the  errors  can  be  reduced  to  the  point  where  the  across- 
shelf  geostrophic  current  variations  are  well  resolved.  Recall 
that,  in  the  arbitrary  case,  the  density  variation  is  dominated  by 
salinity  but  reinforced  by  the  temperature  variation.  This  may 
be  contrasted  with  the  weak  compensation  case  [Fig.  10(a)], 
for  which  the  errors  in  the  computation  of  dD/dx ,  dVg/dz , 
and  Vg(z0)  (Fig.  10(b)-(d),  respectively)  almost  swamp  the 
signal,  in  spite  of  the  spatial  filtering.  Thus,  the  geostrophic 
current  calculation  will  succeed  only  in  certain  regions  and 
seasons  when  the  SSD  signal-to-noise  ratio  is  sufficiently  high 
to  resolve  the  across-shelf  variability. 


Fig.  10.  Derivation  of  geostrophic  current,  as  in  Fig.  9,  but  for  weakly 
compensating  SST  and  SSS  variation  (the  other  parameters  are  as  in  Fig.  9). 
The  transects  shown  are  for  surface  (a)  density,  (b)  horizontal  density  gradient, 
and  (c)  geostrophic  shear,  and  (d)  geostrophic  velocity  at  depth  ( zo ),  relative  to 
the  surface.  This  corresponds  to  the  weakly  compensated  transect  of  Fig.  6. 


Fig.  11.  (a)  STARRS  aircraft  with  the  instrument  mounted  underneath, 

(b)  instrument  system  comprising  L-,  C-,  and  IR-band  radiometers,  and 

(c)  L-band  instrument  six-beam  sampling  pattern.  The  average  pixel  size  is 
1  km  at  a  nominal  2600-m  flight  altitude. 

V.  Applications 

STARRS  [Fig.  11]  is  an  airborne  radiometer  system  de¬ 
signed  to  map  SSS,  SST,  and  SSR  with  a  spatial  resolution  of 
about  1  km  over  rectangular  domains  with  linear  dimensions 
of  about  60  km  in  a  4-h  flight.  (See  Appendix  A-l  for  a 
description  of  the  instrument  and  its  calibration,  and  [30]  for 
details  on  instrument  performance.)  Surveys  of  the  Mississippi 
River  plume  were  conducted  using  STARRS  in  the  early  boreal 
summer  in  May  and  in  the  early  winter  in  November  2004 
[37].  Surveys  of  the  Plata  River  plume  were  conducted  in  the 
austral  winter  in  August  and  September  2003  and  in  the  summer 
between  January  30  and  February  10,  2004.  (See  [38],  [39], 
and  Appendix  A-2  for  details  on  the  survey  operations  and  data 
processing.)  Only  SSS  and  SST  were  successfully  surveyed, 
as  the  C-band  instrument  designed  to  measure  SSR  suffered 
excessive  calibration  noise  (since  remedied),  so  the  roughness 
effects  could  not  be  directly  observed.  Instead,  wind  speed 
observations  from  meteorological  buoys  and  ships  were  used  to 
perform  roughness  corrections  based  on  the  Camps  et  al.  [33] 
empirical  formula.  The  survey  data  were  used  here  to  study  the 
seasonal  variation  in  density  distribution,  as  illustrated  by  the 
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Fig.  12.  Seasonal  density  contrasts  for  the  Mississippi  River  plume  from  STARRS  SST,  SSS,  and  SSD  maps  for  (a),  (c),  and  (e)  winter  and  (b),  (d),  and  (0 
summer  with  T s  and  S  effects  compensating  and  reinforcing,  respectively.  In  both  cases,  S  dominates  T  in  its  effects  on  D.  The  box  encompasses  the  area 
occupied  by  the  SEED  mooring  array.  Asterisks  mark  the  respective  northwest  and  southeast  locations  of  NOAA  data  buoys  42  007  and  42  040.  The  IR  radiometer 
nadir  view  of  SST  (a)  is  repeated  across  track,  and  all  swath  widths  are  exaggerated  by  a  factor  of  2  for  clarity. 


results  obtained  for  the  Mississippi  plume,  and  the  across-shelf 
variation  of  near-surface  geostrophic  velocity,  as  illustrated 
using  the  Plata  plume  results. 


A.  Seasonal  Density  Contrasts 

Representative  STARRS  SST  survey  results  for  the 
Mississippi  River  from  the  winter  flight  on  November  7,  2004 
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Fig.  13.  STARRS  survey  spanning  the  continental  shelf  off  Brazil  that  was  flown  in  winter  2003,  showing  a  portion  of  the  Plata  plume  evident  as  a  band  of 
coastal  water  low  in  (a)  SST,  (b)  SSS,  and  (c)  SSD.  (d)  A  ship  survey  conducted  four  days  earlier  along  the  Albardao  showing  a  vertical  section  of  D  (ere). 


[Fig.  12(a)]  showed  cool  waters  (~23  °C)  near  shore  in  the 
vicinity  of  the  Mississippi  Delta  (29°  S,89°  W)  and  at  the 
western  entrance  into  Mississippi  Sound  (30°  S,89°  W)  and 
warmer  waters  (~27  °C)  seaward  of  these  locations.  This 
may  be  compared  with  the  earlier  summer  flight  on  May  7, 
2004.  The  summer  SST  survey  [Fig.  12(b)]  shows  warmer 
temperatures  (~26  °C)  pervading  the  coastal  domains  of  the 
delta  and  entrance  to  Mississippi  Sound  with  a  seaward  exten¬ 
sion  into  the  western  Mississippi  Bight.  These  contrast  with 
cooler  offshore  temperatures  (~24  °C)  south  and  east  of  the 
delta. 

The  winter  and  summer  SSS  maps  [Fig.  12(c)  and  (d)]  show 
similar  patterns,  with  fresher  waters  near  shore  and  more  saline 
water  well  south  and  east  of  the  delta.  While  the  extreme 
salinity  values  are  similar  (5  and  30  psu),  the  area  of  freshwater 
influence  is  more  extensive  in  the  summer.  However,  during 
that  season,  a  streamer  of  fresher  water  moved  eastward  and 
southward  from  the  delta  and,  in  subsequent  surveys  (not 
shown),  entered  the  western  end  of  NRL’s  “Slope  to  Shelf  Ener¬ 
getics  and  Exchange  Dynamics  (SEED)”  project  mooring  array, 
as  confirmed  by  shipboard  CTD  profiles  [37].  The  salinities 
observed  by  STARRS  agree  with  near-surface  (~2-m  depth) 
CTD  measurements  obtained  inside  the  array,  and  the  surface 


temperatures  are  consistent  with  those  from  NOAA  National 
Data  Buoy  Center  data  buoy  42  040  in  the  same  area. 

The  winter  and  summer  SSD  distributions  [Fig.  12(e)  and 
(01.  respectively,  closely  follow  the  corresponding  SSS  dis¬ 
tributions,  which  dominate  the  density  variability,  i.e.,  lower 
density  water  appears  around  the  delta  and  in  the  entrance  to  the 
sound,  with  a  seaward  extension  into  Mississippi  bight  evident 
in  the  summer.  However,  the  influence  of  the  contrasting  winter 
and  summer  SST  distributions  is  evident,  particularly  around 
the  delta  and  near  the  coast  to  the  north,  with  SSD  showing 
a  weaker  density  gradient  in  these  areas  in  the  winter.  This 
suggests  that  reinforcement  in  summer  has  intensified  the  den¬ 
sity  contrast  in  comparison  with  the  compensating  situation  in 
winter.  It  is  notable  that,  in  contrast  to  winter,  the  density  is 
lower  within  the  mooring  array  in  the  summer  despite  SST’s 
being  cooler  locally.  This  is  associated  with  the  lower  salinities 
and  is  a  further  indication  of  salinity  influence  dominating 
temperature  throughout  the  domain. 

B.  Density- Induced  Geostrophic  Current 

The  wintertime  STARRS  survey  (Legs  I— III)  conducted  off 
the  Atlantic  coast  of  South  America  showed  the  Plata  plume 
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Fig.  14.  STARRS  data  from  the  Albardao  transect  (shown  in  Fig.  13)  plotted 
as  an  across-shelf  profile  through  the  Plata  plume  showing  (a)  SST,  (b)  SSS, 
(c)  SSD,  (d)  horizontal  density  gradient,  (e)  Turner  vector,  and  (f)  geostrophic 
velocity  at  the  base  of  a  50-m-deep  surface  layer. 


extending  from  the  Rio  de  la  Plata  entrance  (latitude:  36°  S) 
to  north  of  Florianopolis  (28°  S)  along  the  Uruguayan  and 
southern  Brazilian  continental  shelf  [38],  [39].  Leg  II  of  the 
survey  showed  cool  fresh  low-density  water  [Fig  1 3(a) — (c)] 
spread  along  the  inner  shelf,  with  warmer,  more  saline,  and 
denser  water  seaward.  The  density  is  low  (~17  at)  along 
the  inner  continental  shelf  and  higher  in  the  outer  shelf  and 
slope,  increasing  to  ~23  at  in  the  northern  part  of  the  domain 
and  to  ~25  at  in  the  south,  with  a  maximum  across-shelf 
contrast  of  8  ct  units.  This  across-shelf  SSD  contrast  is  most 
readily  seen  in  plots  along  the  Albardao  transect  [Fig.  14]. 
As  demonstrated  for  the  Mississippi  plume  in  the  previous 
section,  this  high-density  contrast  resulted  from  the  dominating 
effect  of  salinity  due  to  the  Plata  river  plume  [Fig.  14(b)] 
acting  in  spite  of,  and  partially  compensated  by,  the  generally 
cooler  waters  near  the  coast  [Fig.  14(a)].  The  salinities  and 
densities  are  similar  to  those  measured  by  the  ship  on  the  outer 
shelf  four  days  earlier  but  are  lower  near  shore.  The  STARRS 
temperatures  are  higher  by  an  average  of  2.5  °C,  and  the 
coastal  salinities  are  lower  and/or  extend  further  seaward.  These 
differences  could  be  explained  by  the  subsequent  ship  weather 
observations  obtained  admittedly  5°  of  latitude  (~500  km) 
further  north,  which  show  that  the  air  temperature  warmed  by 
3  °C  while  the  humidity  rose  by  10%  and  the  winds  shifting 
from  south-southwest  to  west  became  less  favorable  for  down- 
welling.  In  addition,  STARRS  measures  SSS  in  the  top  few  cen¬ 
timeters  of  the  ocean,  whereas  the  ship  SSS  measurements  are 
from  the  bulk  (~2-m  depth).  Thus,  STARRS  will  likely  observe 
fresher  salinities  in  the  presence  of  a  thin  or  stratified  buoyant 
plume. 

A  dimensional  analysis  of  the  along-  and  across-shelf  mo¬ 
mentum  equations  [38]  showed  that  the  Plata  river  plume,  in 
winter,  is  approximately  in  semigeostrophic  balance,  with  wind 
dominating  bottom  drag  in  the  along-shelf  (ageostrophic)  bal¬ 


ance  and  across-shelf  pressure  gradient  balancing  Coriolis  force 
in  the  across-shelf  (geostrophic)  balance.  We  use  this  result 
to  justify  using  the  “thermal  wind”  term  in  (1)  to  investigate 
the  effects  of  geostrophic  shear  in  the  near-surface  waters  over 
the  shelf.  The  shipboard  density  vertical  section  obtained  along 
the  Albardao  transect  (Fig.  13(d),  plotted  from  [40])  confirms 
the  STARRS  SSD  observation  of  increasing  density  seaward 
across  the  shelf  and  particularly  along  the  Albardao  transect. 
This  density  section  also  demonstrates  that  the  across-shelf  den¬ 
sity  gradient  is  relatively  uniform  as  a  function  of  depth,  down 
to  a  level  of  about  50  m,  which  corresponds  to  the  inner  portion 
of  the  continental  shelf  domain.  The  density  sections  along 
the  Rio  Grande  and  Solidao  transects  lying  further  north  (not 
shown)  exhibited  a  similar  uniform  horizontal  density  gradient 
and  depth  dependence.  This  structure  likely  exhibited  weaker 
isopycnal  tilts  when  the  STARRS  survey  was  undertaken  four 
days  later,  in  response  to  wind  shifting  to  the  west,  but  the 
relatively  uniform  pattern  of  horizontal  and  vertical  density 
variations  observed  by  the  ship  is  likely  to  be  a  robust  feature 
of  the  wintertime  Plata  plume.  The  corresponding  STARRS 
horizontal  SSD  profile  along  the  Albardao  section  [Fig.  14(c)], 
which  was  computed  from  the  observed  SST  [Fig.  14(a)]  and 
SSS  [Fig.  14(b)],  shows  a  spatially  fluctuating,  but  generally 
increasing,  density  offshore  (with  the  trend  being  marked  by 
the  solid  straight  line).  The  density  gradient  across  the  plume 
[Fig.  14(d)]  was  mostly  compensating  in  T  and  S  with  S 
variations  dominating  the  weak  T  gradient,  except  where  the 
S  gradient  leveled  off  [Fig.  14(e)]. 

Based  on  the  assumptions  (previously  justified)  that  the 
across-shelf  momentum  balance  is  geostrophic  and  that  the 
horizontal  density  gradient  is  uniform  as  a  function  of  depth 
within  the  specified  depth  range,  we  use  the  geostrophic  theory 
in  Section  IV  and  Appendix  C  to  determine  the  corresponding 
geostrophic  velocity  relative  to  the  surface,  at  the  base  of 
this  50-m-deep  layer.  The  along-shelf  geostrophic  velocity 
[Fig.  14],  which  was  computed  from  the  horizontal  density  gra¬ 
dient  [Fig.  14(d)],  after  low-pass  filtering  density  [Fig.  14(c)] 
using  a  Hanning  filter  with  an  effective  cutoff  wavelength  of 
30  km,  fluctuates  within  the  range  of  —1. 0-0.5  ms-1.  The 
geostrophic  current  at  the  base  of  the  near-surface  layer  was 
strongest  (southward),  where  the  horizontal  density  gradient 
was  dominated  by  salinity.  It  was  reduced  by  the  effects  of 
density  compensation  in  the  middle  and  outer  continental  shelf 
(where  the  Turner  angle  approaches  90°).  It  was  reversed 
where  reinforcement  produced  a  negative  density  gradient 
[Fig.  14(e)],  but  the  contributions  of  the  T  and  S  variations  to 
the  density  gradient  are  weak  in  this  area,  as  evidenced  by  the 
small  Turner  radius,  and  the  actual  density  gradient  is  subject 
to  greater  uncertainty. 

Filtering  was  performed  to  minimize  the  potential  impact  of 
small  spatial  scale  fluctuations  due  to  ageostrophic  processes 
acting  within  the  mixed  layer  over  space  and  time  (inertial  or 
finer)  scales  that  are  too  small  to  be  significantly  affected  by 
the  Earth’s  rotation.  The  filter  cutoff  wavelength  was  chosen 
based  on  estimates  of  the  baroclinic  Rossby  radius  Lr  by 
Chelton  et  al.  [41].  They  mapped  Lr  globally  using  clima¬ 
tological  T  and  S  data  in  1°  x  1°  cells,  with  the  result  that 
Lr  ~  30  km  off  the  South  American  Atlantic  coast.  Given  the 
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complexity  of  the  circulation  in  this  region  and  the  strong  influ¬ 
ence  of  river  runoff  on  stratification,  higher  resolution  estimates 
from  representative  current  meter  moorings  or  high-resolution 
shipboard  acoustic  Doppler  and  hydrographic  measurements 
are  needed  to  determine  the  space  and  time  scales  at  which 
the  geostrophic  theory  is  valid  over  the  continental  shelf.  As 
a  more  conservative  estimate,  which  assumes  a  representative 
spatial  scale  corresponding  to  the  width  of  the  continental 
shelf,  we  also  computed  the  spatial  mean  geostrophic  velocity 
(based  on  the  mean  horizontal  density  gradient  represented  by 
the  horizontal  line  in  Fig.  14(d)).  We  found  a  mean  value  of 
—22  cm  •  s_1,  with  the  minus  sign  indicating  southward  flow 
at  a  depth  of  50  m,  relative  to  the  surface.  In  doing  this,  we 
recognize  that  significant  information  on  geostrophic  current 
structure  could  reside  at  smaller  scales. 

As  mentioned  in  Section  IV,  given  an  estimate  of  the 
barotropic  surface  current  relative  to  a  mean  sea  surface  or  a 
geoid,  the  baroclinic  subsurface  current  computed  here  could 
be  added  to  determine  the  total  geostrophic  current  throughout 
the  near-surface  layer.  Estimation  of  the  barotropic  compo¬ 
nent  of  the  surface  geostrophic  current  using  these  means  has 
already  been  demonstrated  by  others  [  14]— [16]  and  [42]  and 
lies  outside  our  present  scope.  However,  Strub  et  al.  [15] 
discuss  resolution  scales  for  surface  geostrophic  currents  ob¬ 
tained  from  altimetry,  and  Powell  and  Leben  [43]  have  pro¬ 
posed  a  filter  that  optimizes  the  tradeoff  between  reduc¬ 
ing  altimeter  noise  and  maintaining  high  spatial  resolution. 
Strub  et  al.' s  results  suggest  that  surface  geostrophic  velocity 
Vs  can  be  estimated  from  TOPEX  altimeter  data  with  a  preci¬ 
sion  (rms  difference)  of  about  7—8  cm  •  s"1  using  a  50-km-long 
“loess”  filter  and  a  gradient  operator  that  spans  62  km.  With 
this  processing  configuration,  the  cross-track  velocities  were 
found  to  resolve  features  with  horizontal  scales  in  the  range 
of  50-80  km.  The  along-track  spatial  resolutions  of  satellite- 
borne  L-band  salinity  mapping  instruments  (SMOS  and 
Aquarius)  over  the  deep  ocean  are  expected  to  be  comparable 
to  those  of  the  altimetric  data  [44].  This  suggests  the  possibility 
of  combining  data  from  the  altimeter  and  L-band  radiometer 
systems  to  estimate  absolute  geostrophic  currents  near  sur¬ 
face  in  deep  ocean  regions  at  spatial  scales  on  the  order  of 
100  km.  Another  possibility  is  suggested  by  a  recent  study  by 
Isem-Fontanet  et  al.  [17],  which  employed  microwave  SST 
data  of  ~50-km  resolution  to  determine  density-driven  near¬ 
surface  currents  using  an  approximate  SSS  value  and  the 
QG  theory.  This  theory  is  likely  to  be  valid  where  mixed 
layer  deepening  or  strong  wind  events,  as  in  winter,  produce 
a  close  relationship  between  surface  density  anomalies  and 
the  circulation  in  which  they  are  embedded.  The  derivation 
of  SSD  from  L-band  microwave  radiometry,  as  demonstrated 
here,  could  complement  the  QG  theory  and  expand  its  do¬ 
main  of  application  to  situations  where  knowledge  of  the  SST 
alone  is  not  sufficient  to  determine  the  SSD  anomaly  field. 
The  QG  method  and  our  thermal- wind-based  method  share 
an  advantage  over  the  barotropic  methods  in  that  they  relate 
the  horizontal  density  structure  to  the  vertical  velocity  struc¬ 
ture  and  are  capable  of  providing  estimates  of  the  vertical 
shear,  and  also  the  internal  velocities,  if  the  surface  velocity  is 
known. 


VI.  Discussion 

We  have  presented  a  simple  algorithm  for  determining  den¬ 
sity  from  the  microwave  brightness  and  physical  temperature  of 
the  sea  surface  based  on  the  nonlinear  Klein  and  Swift  emissiv- 
ity  relation  and  the  UNESCO  SES  [7],  [11].  The  nonlinearities 
result  in  a  sensitivity  of  errors  in  D  to  errors  in  the  observed 
Ts  and  Tb  that  significantly  varies  throughout  the  Ts-Tb 
space.  These  errors  were  quantified  for  several  idealized  cases 
in  which  the  effects  of  T  and  S  either  mutually  compensate 
or  reinforce  one  another  in  their  effect  on  S.  The  analysis  was 
extended  to  determine  the  expected  precision  and  bias  of  SSD 
retrievals  from  the  STARRS  airborne  radiometer  system  and 
from  the  SMOS  and  Aquarius  salinity  mapping  satellites.  We 
then  applied  a  density  compensation  index,  the  “Turner  vector” 
based  on  the  density  ratio  [  1 8],  and  the  horizontal  Turner  angle 
[19]  to  give  a  compact  representation  of  both  the  significance 
and  influence  of  S  and  T  on  D  and  showed  that  errors  in 
retrieving  horizontal  density  gradients  depended  on  whether 
the  S  and  T  contributions  were  mutually  compensating  or 
reinforcing. 

We  used  the  SSD  algorithm  to  retrieve  SSD  information 
from  SST  and  SSS  observations  obtained  from  the  STARRS 
airborne  IR  and  microwave  radiometers,  with  sufficient  pre¬ 
cision  for  several  practical  oceanographic  applications.  These 
included  mapping  the  surface  density  of  two  major  river  plumes 
spanning  the  continental  shelf,  determining  the  relative  com¬ 
pensation  or  reinforcement  effects  of  S  and  T  on  D,  in  the 
Mississippi  plume  during  summer  and  winter,  and  estimating 
surface  geostrophic  shear  and  subsurface  geostrophic  current  in 
the  Plata  plume  during  winter. 

The  Mississippi  and  Plata  river  plumes  showed  evidence  of 
density  compensation  and  density  reinforcement,  associated 
with  cool  or  warm  seasons  having  offshore  temperature  gradi¬ 
ents  of  the  opposite  sign.  For  the  Mississippi  River  plume,  the 
SST  and  SSS  were  either  correlated  (winter)  or  anticorrelated 
(summer)  and  alternately  showed  compensation  or  reinforce¬ 
ment,  respectively.  In  contrast  to  the  Mississippi  plume,  the 
Plata  plume  was  of  larger  scale,  and  in  winter,  it  formed  a 
persistent  buoyant  coastal  current  that  propagated  northward, 
under  the  combined  influence  of  Earth’s  rotation  and  southerly 
winds. 

The  density  gradient  across  the  Plata  plume  was  mostly 
compensating  in  T  and  S  with  S  variations  dominating  the 
weak  T  gradient,  except  where  the  S  gradient  leveled  off. 
Since,  under  these  conditions,  the  buoyant  plume  exhibited  a 
semigeostrophic  momentum  balance,  it  was  possible  to  use 
the  across-shelf  horizontal  density  gradient  to  compute  the 
geostrophic  current  shear  at  the  surface.  Subject  to  appropriate 
assumptions,  the  subsurface  geostrophic  current  due  to  the  hor¬ 
izontal  density  gradient  was  inferred.  The  assumptions,  which 
allow  for  an  arbitrary  variation  of  density  in  the  horizontal 
direction  and  linear  variation  in  the  vertical  direction  are  not 
overly  restrictive,  and  the  linear  constraint  could  be  relaxed 
when  more  information  is  available  on  the  vertical  density 
variation.  The  current  was  strongest  (southward)  where  the 
horizontal  density  gradient  was  dominated  by  salinity.  It  was 
reduced  by  the  effects  of  density  compensation  in  the  middle 
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and  outer  continental  shelf  and  reversed  where  reinforcement 
produced  a  negative  density  gradient. 

VII.  Conclusion 

Based  on  the  foregoing  algorithm  development,  assessments 
of  algorithm  performance,  and  our  experience  with  STARRS, 
the  determination  of  SSD  from  satellite  IR  and  microwave 
observations  of  SST  and  SSS  appears  feasible.  More  work 
is  needed  to  identify  the  likely  error  levels  (both  accuracy 
and  precision)  for  SSD  determined  from  the  satellite  salinity 
mapping  technology,  under  conditions  that  will  be  encountered 
at  the  surface  of  the  deep  ocean.  The  estimation  of  geostrophic 
currents  associated  with  horizontal  density  gradients,  which 
could  be  added  to  altimetric  surface  geostrophic  current  veloc¬ 
ity  estimates  to  determine  the  near-surface  geostrophic  current 
profile,  can  also  be  envisaged  under  appropriate  assumptions 
and  in  specific  geographic  regions. 

This  paper  has  focused  on  the  retrieval  of  SSD  and 
geostrophic  shear  from  airborne  radiometer  systems  in  coastal 
regions  subject  to  buoyancy-driven  coastal  currents.  Future 
studies  will  focus  on  determining  the  feasibility  of  the  SSD 
and  geostrophic  shear  retrieval  from  the  SMOS  and  Aquarius 
satellite  systems  in  specific  regions  and  seasons  subject  to  mod¬ 
erate  to  strong  buoyancy  contrasts,  such  as  western  boundary 
currents,  and,  more  generally,  at  the  margins  of  the  synoptic 
and  basin  scale  where  the  signals  are  likely  weaker  but  broader 
in  the  horizontal  scale. 

Our  error  estimates  suggest  that  SSD  can  be  determined  in 
temperate  and  tropical  climates  at  approximately  equal  levels 
of  precision  and  accuracy  of  ~0.75  and  0.25  kg  •  m~3  for 
SMOS  and  Aquarius,  respectively,  with  temporal  but  no  spa¬ 
tial  averaging,  at  resolutions  of  ~  100  km.  If  the  postlaunch 
error  levels  for  satellite  determinations  of  both  surface  density 
SSD  and  vertical  shear  of  the  geostrophic  velocity  dVg/dz 
prove  to  be  acceptable,  density  products  derived  from  these 
satellites  could  be  developed  and  used  in  a  number  of  ways. 
For  example,  the  SSD  information  could  be  directly  assimi¬ 
lated  into  numerical  hydrodynamic  models  as  an  alternative 
to  the  current  practice  of  independently  specifying  surface 
temperature  and  salinity.  This  might  have  computational  ad¬ 
vantages  for  certain  types  of  models,  specifically  those  that 
were  set  up  to  investigate  dynamical  issues.  Identical  twin 
experiments  could  be  performed  to  determine  the  relative  ad¬ 
vantages  of  assimilating  either  SST  and  SSS  or  SSD  in  such 
models.  The  products  could  be  used  independently  of  the 
models,  or  they  could  exhibit  a  weak  dependence  on  them; 
models  and  in  situ  mooring  data  could  be  used  to  determine 
the  validity  of  the  assumptions  of  geostrophy  and  uniformity 
of  subsurface  horizontal  density  gradient  in  specific  regions. 
In  either  case,  the  production  of  an  operational  algorithm 
that  generates  estimates  of  both  SSD  and  dVg(z)/dz  could 
be  useful  in  combination  with  estimates  of  SSH  and  sur¬ 
face  geostrophic  current  Ve  available  from  satellite  altime¬ 
ters.  This  suggests  that  existing  products  such  as  OSCAR 
[16],  [45],  which  account  for  both  geostrophic  and  Ekman 
effects  to  estimate  surface  drift,  could  be  upgraded  to  take 
advantage  of  the  baroclinic  shear  deduced  from  the  salinity 


mapping  satellites  to  project  surface  current  estimates  into  the 
subsurface. 

Appendix  A 

Instrumentation  and  Data  Processing 

I)  Instrumentation :  STARRS  (Fig.  11)  is  an  airborne  ra¬ 
diometer  system  built  for  NRL  to  map  SSS  in  estuaries  and 
marginal  seas.  It  is  flown  beneath  a  twin-engined  aircraft 
[Fig.  1 1(a)]  and  comprises  L-  and  C-band  microwave  radiome¬ 
ters  and  an  IR  radiometer  [Fig.  11(b)]  with  an  integrated  Global 
Positioning  System  (GPS)  receiver  and  fiber-optic  gyro  [30], 
[39],  [46].  Its  primary  purpose  is  to  map  coastal  SSS,  SST,  and 
SSR.  The  L-band  radiometer,  which  is  used  to  map  SSS,  is  a 
multibeam  system  sensing  natural  microwave  emission  from 
the  sea  surface  within  a  24-MHz-wide  protected  band  centered 
at  a  frequency  of  1.413  GHz.  It  has  a  brightness  temperature 
sensitivity  of  0.35  K  at  a  sampling  interval  of  2  s.  This  is 
equivalent  to  0.52  K  at  1  s  and  can  be  reduced  to  0.28  and 
0.21  K  by  averaging  over  1 2  and  24  s,  respectively.  Due  to  some 
sample  autocorrelation,  this  reduction  is  less  than  Gaussian 
(which  scales  as  \/y/~N>  where  N  is  the  number  of  samples 
averaged).  The  six  antenna  beams,  each  with  15°  half-power 
beam  width,  point  downward  and  to  either  side  of  the  aircraft 
at  incidence  angles  of  ±7°,  22°,  and  38°  [Fig.  11(c)].  The 
C-band  radiometer,  which  can  sense  both  SSR  and  SST,  has 
a  single  20°  wide  beam  directly  pointing  beneath  the  aircraft 
(nadir  viewing).  It  senses  natural  emission  within  six  channels, 
at  5.2,  5.6,  5.9,  6.2,  6.6,  and  7.1  GHz.  The  nadir-viewing  IR 
radiometer  measures  SST  independently  from  thermal  emission 
in  the  8-14-/xm  and  9.6-1 1.5-/im  bands.  In  contrast  to  the 
L-  and  C-band  radiometers,  it  cannot  penetrate  clouds,  so 
gaps  in  the  IR  observations  are  filled  by  interpolation.  The 
GPS  and  gyro  provide  Coordinated  Universal  Time  (UTC) 
time  and  position,  and  aircraft  attitude,  respectively.  For  typ¬ 
ical  aircraft  altitudes  of  2600  m,  the  six  L-band  beams  have 
footprints  with  a  width  of  700-1100  m  and  a  swath  of 
5.2-km. 

STARRS  instruments  are  calibrated  before  and  after  each 
mission.  The  IR  radiometer  is  calibrated  at  various  steps  span¬ 
ning  the  observational  temperature  range  using  a  blackbody 
target  immersed  in  a  temperature  controlled  water  bath.  The 
L-  and  C-band  instruments  are  more  difficult  to  calibrate 
due  to  their  large  antenna  size  and  wide  Tb  range.  This  is 
done  first  in  the  laboratory  by  observing  a  room  temperature 
(~300  K)  target  comprising  a  microwave  absorbing  mat,  with 
its  temperature  monitored  to  provide  a  “warm"  calibration 
point.  The  instrument  is  then  set  up  outside  at  night  (avoiding 
solar  radiation)  to  observe  the  brightness  temperature  of  the 
cold  sky  (~8  K),  which  provides  the  “cold"  calibration  point. 
The  sky  temperature  is  determined  by  the  constant  cosmic 
background  plus  galactic  emissions  in  the  field  of  view  and 
is  augmented  by  atmospheric  O2  emission  and  attenuated  by 
atmospheric  absorption.  Instrument  stability  is  assured  by  ther¬ 
mal  regulation  and  the  measurement  of  internal  temperatures 
of  key  components,  including  the  antennas,  and  hot  and  cold 
calibration  reference  targets.  Calibration  involves  regressing  the 
observed  brightness  temperature  of  the  external  and  internal 
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targets  against  the  known  target  brightness  temperatures  while 
taking  into  account  the  internal  physical  temperatures  of  the 
components  and  reference  noise  sources  to  obtain  a  multivariate 
calibration  curve.  Internal  sampling  and  processing  operations 
minimize  the  effects  of  instrument  drift  and  gain  variations  and 
temporally  average  the  data  to  maximize  sensitivity  [30].  Dur¬ 
ing  operational  deployment,  the  internal  calibration  targets  are 
sampled  at  regular  intervals  and  used  to  correct  the  interspersed 
Tb  observations  of  the  seawater  surface.  These  vary  over 
~  100-130  K,  depending  on  the  incidence  angle  and  allowing 
for  the  various  factors  mentioned  in  Section  II. 

2)  Data  Acquisition  and  Processing:  Surveys  of  the 
Mississippi  River  plume  using  STARRS  were  conducted  in 
early  boreal  summer  during  May  2-8  and  in  early  winter 
during  November  4-17,  2004.  The  flights  coincided  with  in 
situ  shipboard  sampling  of  T  and  5,  and  were  bracketed  by 
the  deployment  of  an  extensive  mooring  array  spanning  the 
continental  shelf  break  east  of  the  Mississippi  Delta,  as  part 
of  NRL’s  SEED  experiment.  Surveys  of  the  Plata  River  plume 
were  conducted  in  the  austral  winter  between  August  19  and 
September  5,  2003  and  in  summer  between  January  30  and 
February  10,  2004.  The  flights  spanned  the  continental  shelf 
and  extended  from  northern  Argentina  (latitude:  38°  S)  to 
Southern  Brazil  (latitude:  27°  S),  which  is  a  distance  of  some 
1 200  km,  including  the  outflow  from  the  Rio  de  la  Plata  estuary 
and  Patos  Lagoon  (see  [39]  for  a  more  detailed  description  of 
the  mission  and  salinity  survey  results).  The  flights  reported 
here  were  performed  in  the  line-transect  mode  with  tracks 
separated  by  gaps,  in  contrast  to  the  mapping  mode,  which  gives 
complete  spatial  coverage.  Altitudes  were  chosen  to  optimize 
the  sampling  subject  to  operational  constraints,  such  as  aircraft 
range  (duration)  and  cloud  cover. 

The  calibrated  STARRS  thermal  IR  radiometer  data  were 
used  to  determine  SST  by  using  a  dual-window  algorithm 
with  corrections  for  the  effects  of  the  intervening  water  va¬ 
por  obtained  from  estimates  of  precipitable  water  vapor.  The 
latter  water  vapor  was  estimated  using  surface  measurements 
of  wet  and  dry  bulb  temperature  (or  dew  point)  from  me¬ 
teorological  buoys  (during  the  SEED  surveys)  or  the  ship’s 
weather  station  (during  Plata),  and  an  assumed  exponential 
falloff  in  water  vapor  concentration  with  height  [47],  with 
a  22-km  scale  height  [48,  p.  171].  The  calibrated  L-band 
brightness  temperatures  were  first  edited  to  remove  anomalous 
pixels  due  to  the  presence  of  land  in  the  field  of  view  or 
aircraft  maneuvers,  and  data  obtained  from  beam  3R  during  the 
May  2004  survey  were  omitted  due  to  an  instrument  noise 
problem  affecting  only  that  beam,  which  has  since  been  reme¬ 
died.  The  remaining  Tb  values  were  corrected  for  reflected 
atmospheric  02  and  extraterrestrial  (cosmic  and  galactic)  emis¬ 
sion,  atmospheric  absorption,  and  wind-induced  roughness  and 
then  converted  to  SSS  by  empirical  inversion  of  the  KS77 
relation.  The  roughness  correction  was  performed  using  the 
Camps  et  al.  [33]  empirical  formula  using  representative  wind 
speeds  from  buoy  and/or  ship  anemometers  adjusted  to  the 
standard  meteorological  height  (10  m).  The  KS77  algorithm  ac¬ 
counted  for  the  effects  of  SST,  polarization,  and  beam  incidence 
angle  (coupled  with  aircraft  pitch  and  roll)  on  the  retrieved 
SSS.  An  additional  empirical  roll  bias  correction  was  applied  to 


correct  across-track  brightness  temperature  differences:  1.5°  in 
the  November  2004  data  (SEED)  and  0.5°  in  the  August  (Plata) 
data.  These  differences  were  likely  due  to  gyro  initialization 
errors  since  subsequent  gyro  replacement  eliminated  this  error. 
Residual  across-track  differences  evident  in  the  November  2004 
data  are  likely  due  to  small  beam-to-beam  calibration  biases 
that  are  consistent  with  the  known  instrument  Tb  calibration 
precision  (  0.7  K). 


Appendix  B 
Turner  Vector 


The  density  ratio  defines  the  relative  contribution  of  thermal 
and  haline  gradients  and  associated  anomalies  to  the  observed 
density  gradients.  It  is  defined  as  [18] 


t  _  adT/dx 
p  “  pdS/dx 


(3) 


where  a  and  are  the  coefficients  of  thermal  expansion 
and  haline  contraction,  respectively.  These  are  defined  as 
[13,  p.  254] 


In  (3),  we  express  the  anomalies  as  spatial  gradients  in  T  and 
S  as  a  function  of  horizontal  distance  x  offshore,  in  contrast  to 
the  usual  representation  of  Rp  used  in  fine-structure  studies  (in 
which  gravity  plays  a  key  role),  in  terms  of  vertical  derivatives 
in  depth  2. 

Using  the  definition  of  the  total  derivative  (on  the  left) 
with  respect  to  partial  derivatives  (on  the  right),  the  horizontal 
density  gradient  is 


dD/dx  = 


( dD\  (dT\  fdD\  fdS\ 
\dr)TP\dx)  +  \ds )TP\dx) ' 


(5) 


This  can  be  expressed  in  terms  of  a  nondimensional  density 
and  the  preceding  quantities  (4)  as 


1  dD  _  dT  _  dS 
D  dx  a  dx  ^  dx 


(6) 


If  S  and  T  variations  are  both  significant  (in  the  sense  de¬ 
fined  here),  dD/dx  =  0  implies  perfect  compensation,  whereas 
if  dD/dx  is  large,  either  reinforcement  is  occurring  or  one 
of  T  and  S  is  strongly  dominating  D.  To  find  out  whether 
compensation  or  reinforcement  is  occurring  and  which  of  T 
or  S  dominates  the  density  gradient,  we  can  examine  the  sign 
of  the  density  ratio  Rp  and  its  magnitude,  or  absolute  value 
\RP\  relative  to  unity.  If  Rp  >  0,  T  and  S  compensate  in  their 
effect  on  D,  whereas  if  Rp  <  0,  they  reinforce  their  effect  on 
D.  If  \RP\  >  1,  T  effects  dominate  those  of  5,  and  if  \RP\  <  1, 
S  dominates  T.  To  determine  the  utility  of  this  approach,  we 
plotted  Rp  and  \RP\  for  the  arbitrary  density  profile  of  Fig.  4(c) 
and  found  that  reinforcement  occurs  most  of  the  way  across 
the  shelf,  whereas  compensation  occurs  at  the  outer  end.  In 
addition,  S  dominated  T  everywhere,  except  in  a  narrow  band 
at  D  =  40  km,  where  Rp  changed  sign. 
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As  previously  implied,  if  dD/dx  «  0,  this  could  be  due  to 
perfect  compensation  of  T  and  S,  but  this  condition  could  also 
be  met  if  T  and  S  show  no  significant  variation  ( dT/dx  «  0, 
dS/dx  «  0).  In  this  case,  Rp  is  a  ratio  of  two  very  small 
numbers  subject  to  noise  and  is  thus  poorly  defined.  To  reveal 
such  cases,  we  need  an  indicator  of  the  absolute  strength  of 
the  corresponding  variations  in  T  and  S.  An  inconvenient 
singularity  also  arises  in  Rp  when  dS/dx  =  0,  causing  it  to 
become  infinite.  An  approach  that  both  avoids  the  singularities 
and  properly  represents  the  four  combinations  arising  from 
each  pair  of  compensation  and  dominance  states  was  proposed 
by  [18]  in  the  form  of  a  four-quadrant  arctangent  “Turner 
angle”  Tu  expressed  in  degrees,  which  is  related  to  the  density 
ratio  by 


Rp  =  —  tan(Tix  +  45).  (7) 


A  simple  four-quadrant  arctangent  of  the  terms  in  the  density 
ratio  would  serve  our  purposes  for  studies  of  horizontal  den¬ 
sity  variations.  However,  we  chose  to  adopt  the  slightly  more 
esoteric  Turner  angle  (with  its  45°  offset)  for  consistency  with 
studies  of  double  diffusive  processes,  which  focus  on  vertical  T 
and  S  gradients. 

Ruddick  [18]  defines  Tu  in  terms  of  Cartesian  coordinates 
corresponding  to  appropriately  signed  values  of  the  density 
ratio  numerator  and  denominator,  i.e.. 


Ns  = 


(8) 


In  his  notation,  the  signs  are  chosen,  so  positive  (7V|,./Vj.) 
corresponds  to  gravitationally  stable  S  and  T  stratification 
(dp/dz  >  0,  with  2  pointing  down).  For  application  to  horizon¬ 
tal  variations  [19],  for  which  we  substitute  the  x  for  2,  gravity 
is  unimportant,  and  the  signs  merely  determine  the  x-direction, 
which  is  considered  to  represent  a  positive  density  gradient. 
We  choose  the  signs  in  (8)  to  be  positive  and  negative  for  the 
first  and  second  terms,  respectively,  so  that  a  positive  density 
gradient  dp/dx  >  0  corresponds  to  D  and  S  increasing,  and  T 
decreasing,  in  the  x-direction.  Ruddick’s  use  of  the  acceleration 
due  to  gravity  g  and  exponent  2  (which  is  symbolic  and  not 
mathematical)  is  evidently  linked  to  the  definition  of  the  square 
of  the  Brunt-Vaisalla  frequency  in  terms  of  vertical  density 
gradient,  which  does  not  concern  us  here  either. 

The  Turner  angle,  which  is  defined  in  these  terms  using  the 
“atan2”  label  of  the  FORTRAN  and  MATLAB  (The  Mathworks 
Inc.)  syntax  to  denote  the  four-quadrant  arctangent  function,  is 

Tu  =  arctan  2  {Nj*  —  TVj,  N§  -f  N £)  .  (9) 

As  an  indication  of  the  combined  strength  of  the  £  and  T  gra¬ 
dients  influencing  D  (regardless  of  whether  they  are  compen¬ 
sating  or  reinforcing  and  which  dominates),  we  further  define 
the  Turner  “radius.”  The  radius,  with  units  of  square  frequency, 
represents  the  combined  magnitudes  of  the  coordinates  defined 
in  (8),  i.e., 


Tr  =  (10) 


The  resulting  vector  (Tr,  Tu),  which,  adapting  Ruddick’s  ter¬ 
minology,  we  term  the  “Turner  vector”  can  be  plotted  either  in 
point  form  on  a  polar  coordinate  frame  or  in  component  form, 
as  a  function  of  the  space  dimension  in  a  Cartesian  frame  (see 
Figs.  1  and  2). 

In  our  polar  representation  (e.g..  Fig.  8),  which  can  be 
obtained  from  that  of  Ruddick  [18,  Fig.  1],  by  rotating  his  disc 
out  of  plane  about  its  Tu  =  0  axis  ( Rp  =  -1  diagonal).  Turner 
angles  falling  within  the  two  right-side  quadrants  correspond  to 
positive  horizontal  density  gradients,  whereas  those  occupying 
the  left-side  quadrants  represent  negative  density  gradients.  In 
quadrants  1  and  3,  T  dominates  S ,  whereas,  in  quadrants  2  and 
4,  S  dominates  T.  If  these  four  quadrants  of  the  circle  are  each 
subdivided  to  form  eight  45°  sectors  numbered  counterclock¬ 
wise  from  the  zero-degree  line  at  the  right,  then  sectors  1, 4,  5, 
and  8  represent  T  and  S  reinforcement,  and  sectors  2,  3,  6,  and 
7  represent  T  and  S  compensation.  Reinforcement  is  strongest 
in  the  0°  and  180°  directions,  whereas  compensation  is  most 
perfect  in  the  ±90°  directions.  Finally,  the  radial  distance  from 
the  center  of  the  circle  represents  the  combined  strength  of  the 
T  and  S  gradients.  In  studies  of  vertical  gradients,  the  right-  and 
left-hand  quadrants  represent  stable  and  unstable  density  strat¬ 
ification,  respectively,  whereas  sectors  2  and  7  represent  salt 
fingering  and  double  diffusive  activity,  respectively.  However, 
neither  characteristic  is  relevant  when  examining  horizontal 
gradients. 


Appendix  C 

Deriving  Geostrophic  Currents 

Geostrophic  current  is  determined  by  the  following  simpli¬ 
fied  horizontal  momentum  balance  equation  (see  standard  texts 
such  as  [49]  or  [50]  for  derivations).  A  right-handed  Cartesian 
coordinate  system  is  used,  with  the  vertical  2-axis  positive 
above  the  mean  sea  surface.  For  momentum  balance  in  the 
x-axis  direction  (oriented  across-shelf  and  positive  offshore  in 
our  STARRS  applications),  the  acceleration  due  to  the  Coriolis 
force  (at  left)  is  canceled  by  that  generated  by  the  pressure 
gradient  force  (at  right),  i.e., 


fV,= 


1  dp 
D(x)  dx 


GO 


Here,  /  =  2fisin(7/)  is  the  Coriolis  parameter,  which  is 
positive/negative  in  the  Northem/Southem  Hemisphere,  Q  is 
the  Earth’s  rotation  rate  (in  radians  per  second),  77  is  the  latitude, 
and  Vg  is  the  speed  of  the  geostrophic  current.  This  flows  in  the 
y-axis  direction  perpendicular  to  the  horizontal  pressure  gradi¬ 
ent  (with  high  pressure  on  the  right  in  the  northern  hemisphere). 
Following  the  Boussinesq  approximation,  the  density,  which 
modifies  pressure  gradient  dP/dx ,  is  represented  by  a  reference 
density  D0,  which  is  here  chosen  to  be  the  observed  surface 
density  D(x). 

The  pressure  gradient  can  arise  from  external  (sea  level) 
or  internal  (density)  differences.  The  vertical  hydrostatic  ap¬ 
proximation  to  the  vertical  momentum  balance  provides  the 
relationship  between  the  density  and  pressure  needed  to 
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determine  the  geostrophic  current  from  density  observations  in 
the  ocean,  i.e.. 


dP/dz  =  —gD(x,z)  (12) 

where  g  is  the  acceleration  due  to  gravity,  which  is  assumed 
to  be  constant  over  the  depths  and  horizontal  locations  of 
interest  here. 

To  maintain  accuracy  when  computing  P  from  (12),  the 
more  general  position-  and  depth-dependent  form  of  the  density 
function  D(x,z)  is  retained,  consistent  with  the  Boussinesq 
approximation.  Given  the  density  variation  with  depth,  (12)  can 
be  integrated  to  give  pressure  at  any  level  z.  We  fix  2  =  0  at 
the  mean  sea  level  but  define  z  =  £(x)  at  the  actual  surface. 
This  allows  for  small  vertical  excursions,  or  more  significantly 
for  geostrophic  balance,  larger  scale  tilt  about  the  mean,  which 
may  horizontally  vary  along  x ,  i.e., 

£(*) 

P(x,z)  =  P(£(x))  +  g  J  D(x,z)dz  (13) 

Z 

where  P(£(x))  is  the  surface  (atmospheric)  pressure. 

If  we  assume  that  density  varies  arbitrarily  in  x  and  linearly 
in  z ,  we  can  write 


D(x,  z)  =  D(x)  (1  —  z/a)  (14) 


where  a  is  the  scale  depth  for  the  vertical  density  variation, 
which  is  empirically  estimated  here.  Substituting  this  expres¬ 
sion  into  (13)  and  the  resulting  P(x,  z)  relation  into  (11) 
yields  a  form  of  the  “thermal  wind”  equation,  which  is  so 
called  in  meteorology  because  the  temperature  effect  on  density 
dominates  vertical  variations  in  geostrophic  wind,  i.e., 

*(*> 

VM—, -/•>“*■  <l5> 

Z 

Moving  D(x)  outside  the  integral  and  then  integrating  the 
depth-dependent  part  between  limits,  we  get 


Vg(z) 


fD(x)  dx 


D( 


x)  U(x) 


£(*)2 

2  a 


(16) 

Differentiation  of  the  product,  allowing  for  the  x-dependence 
of  surface  height,  leads  to 


V9(z)  = 


fD(x) 


{  £>(*)( 


+ 


(J  /' 

=)]}• 

(17) 


d£(x)  £(x)  d£(x) 

dx 

ax? 


) 


Here,  all  fractions  involving  both  £(x)  and  a  can  be  ne¬ 
glected,  and  £(x)  is  dominated  by  2  for  all  depths  of  practical 
interest.  With  this  and  some  further  simplification,  we  get 


Vg(z) 


dD(x) 

dx 


(18) 


If  a  is  known  with  some  certainty,  the  full  expression  (18) 
can  be  used,  but  substituting  a  representative  surface  density 
of  1015  at  with  an  increase  of  10  at  over  100-m  depth  into 
(14)  gives  a  ~  10  km,  so  for  depths  on  the  order  of  100  m, 
e.g.,  over  the  continental  shelf,  the  term  involving  z/a  can  be 
neglected  with  a  velocity  error  on  the  order  of  1%.  With  this 
simplification,  we  have 


v  (,)  =  _  21  dD(x) 

9K  1  f  dx  fD(x)  dx  ' 


(19) 


In  this  final  expression  for  the  total  geostrophic  velocity  at 
depth  2,  the  first  term  on  the  RHS  can  be  recognized  as  the 
surface  geostrophic  velocity  associated  with  the  sea  surface 
slope,  while  the  second  is  the  internal  velocity  relative  to  the 
surface  associated  with  the  horizontal  density  gradient. 
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